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Abstract 

Hydrodynamics can be consistently formulated on surfaces of arbitrary co-dimension in a background 



> 

cn ' space-time, providing the effective theory describing long-wavelength perturbations of black branes. 

r»^ ■ When the co-dimension is non-zero, the system acquires fluid-elastic properties and constitutes what 
f^^ ' is called a fluid brane. Applying an effective action approach, the most general form of the free energy 

o 

en ' constructed to second order in a derivative expansion. This construction generalizes the Helfrich-Canham 
bending energy for fluid membranes studied in theoretical biology to the case in which the fluid is rotat- 
ing. It is found that stationary fluid brane configurations are characterized by a set of 5 elastic response 
coefficients, 5 hydrodynamic response coefficients and 1 spin response coefficient for co-dimension greater 
than one. Moreover, the elastic degrees of freedom present in the system are coupled to the hydrody- 
namic degrees of freedom. For co-dimension- 1 surfaces we find a 13 independent parameter family of 
stationary fluid branes. It is further shown that elastic and spin corrections to (non)-extremal brane 
effective actions can be accounted for by a multipole expansion of the stress-energy tensor, therefore 
establishing a relation between the different formalisms of Carter, Capovilla-Guven and Vasilic-Vojinovic 
and between gravity and the effective description of stationary fluid branes. Finally, it is shown that the 
Young modulus found in the literature for black branes falls into the class predicted by this approach 
- a relation which is then used to make a proposal for the second order effective action of stationary 
blackfolds and to find the corrected horizon angular velocity of thin black rings. 
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1 Introduction 

Higher-dimensional gravity has been shown to be fruitful as a testing ground for theories of hydrody- 
namics. In recent years many new properties, transport and response coefficients of (charged) fluids 
[1, 2, 3, 4, 5, 6] and superfluids [7, 8, 9, 10, 11] have been uncovered through the study of long- wavelength 
fluctuations of black branes. These fluctuations can be along the worldvolume [12, 13, 14, 15] or boundary 
directions [1, 2] where the fluid lives yielding the usual dynamics of viscous fluid flows, or transverse to 
it, originating instead genuine elastic behavior [5, 16, 6]. 

The blackfold approach [17, 18, 19], being the effective hydrodynamic theory that describes the long- 
wavelength perturbations of black branes, has taught us that the fluid system dual to the gravitational 
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object needs not to live on the boundary of the space-time (fluid/gravity) but can also live on the horizon 
(membrane paradigm) or in an intermediate region [20, 15]. In this regime, black holes and black branes 
should be seen as fluid branes [19, 21]: fluids living on dynamical surfaces of arbitrary co-dimension 
embedded in a background space-time. 

The relation between the toy model of a fluid brane and gravitational physics can be thought of in 
the following way. Classically, black holes share many properties with other products of gravitational 
collapse, namely, they are characterized by a very dense distribution of matter but they differ from other 
stellar objects due to the formation of an event horizon. The crust of stars is composed of matter in a 
high pressure state, as the matter distribution that characterizes the black hole, and hence relativistic 
elasticity can be a good approximation for describing deformations of the star crust. It is then intuitively 
understandable why deformations of the matter source of a black hole can give rise to elastic behavior 
[22, 23, 24, 5]. On the other hand, it is also possible in some cases to perform inhomogeneous fluctuations 
of the event horizon without deforming the source leading to an effective viscous fluid behaviour [12, 13, 
14, 15]. These different possibilities are mapped onto the fluid brane toy model as transverse deformations 
to the surface where the fluid lives (elastic), or to fluctuations on the different thermodynamic quantities 
and fluid velocities that characterize the overall fluid conflguration (fluid) [19, 21]. It is then possible to 
think of a black hole as having source degrees of freedom (elastic modes) and horizon degrees of freedom 
(hydrodynamic modes) which in general interact with each other when the entire system is subject to 
a perturbation. The natural question to ask is then: in which situations can the horizon be perturbed 
without perturbing the source and vice-versa? Focusing on stationary black brane conflgurations, in this 
work we partly answer this question when backreaction effects can be neglected. 

Even though fluid branes have gained a role in gravitational physics, the physical problem of how 
a fluid living on a surface reacts to a deformation of the surface is a much older one dating back to 
the first observations of the strange biconcave shape of a red blood cell in the seventeenth century [25]. 
The models that described the shape of a cell, assuming a constant pressure along the surface, are the 
same models that describe the shape of a soap bubble, working under the principle of minimization of 
the surface area. However, in the 60's, Helfrich and Canham proposed what is now called the Helfrich- 
Canham bending energy for fluid membranes [26, 27], consisting of introducing an additional contribution 
to the free energy of the cellular membrane of the form: 

jr[x'^] =a f dAK"^ , (1.1) 



where K is the mean extrinsic curvature vector, dA is the area element induced on the surface and a 
is the modulus of rigidity. The inclusion of this term was enough to understand the biconcave shape of 
the red blood cell [25], a fact that then led to the study of several properties of this fluid-elastic system 
[28, 29, 30, 31, 32]. The rigidity term (1.1) has also found its role in particle physics when Polyakov 
[33], and independently Kleinert [34], added it to the Dirac-Nambu-Goto action in order to obtain an 
improved effective action for QCD. In this context, corrections of the type (1.1) are known as finite 
thickness corrections to brane actions. 



In this paper, we generalize the Helfrich-Canham bending energy to account for possible pressure 
differences along the surface as well as for the possibility of the fluid being in stationary motion. Taking 
a brane-like stringy perspective, what we will be accomplishing is a systematic derivative correction 
to extremal and non-extremal brane actions, including all the possible correction terms quadratic in 
the extrinsic curvature, in the extrinsic twist potential and of second order in derivatives along the 
surface and obtaining the equations of motion for the resulting system. The first covariant formalism 
for obtaining the correct variations of surface tensors was developed by Carter in connection with his 
model of geodynamic-type (or stiff) strings and branes [35, 36, 37]. A more intuitive method that instead 
exploits the symmetries of the worldvolume surface was developed by Capovilla-Guven [38] and applied 
to different physical systems including lipid vesicles [39, 40, 41]. On the other hand, a general framework 
for deriving the equations of motion of finite thickness probe branes based on a multipole expansion of 
the stress-energy tensor was developed by Vasilic-Vojinovic [42, 43]. Carter pointed out in [36] that the 
dynamics of geodynamic-type objects could be accounted for by adding dipole terms to the stress-energy 
tensor. Using mainly the formalism developed by Capovilla-Guven [38] and Vasilic-Vojinovic [43], it is 
shown that this is indeed the case^. 

The upshot of the developments in hydrodynamics in connection with gravitational physics mentioned 
in the beginning of this section is that many systematic methods for constructing theories of fiuid dy- 
namics have been developed [46, 47, 48, 49, 50, 51]. Inspired by partition function methods [46, 47, 48] 
and the effective action approach [49, 50] we construct the free energy for generic stationary fluid flows 
living on dynamical surfaces to second order in a derivative expansion. We show that neutral fluids 
living in a surface of co-dimension greater than one are characterized by 5 elastic response coefficients, 5 
hydrodynamic response coeflicients and 1 spin response coeflicient. Each of these response coefficients is 
associated with a particular term which is added to the free energy and from our analysis we conclude 
that 4 of the elastic terms are coupled to 4 of the hydrodynamic terms due to geometric constraints. 
Contrary to the case of space-filling stationary fiuid fiows which are described by only 3 hydrodynamic 
response coefficients [46], these fiuids are described by a total of 9 independent response coefficients in 
a curved background. Using the generic form of the stress-energy tensor and dipole moment obtained 
from the free energy we match it with the Young modulus measured for black branes [5, 16] and hence 
propose a second order effective action for blackfolds which is then used to predict the corrected horizon 
angular velocity of thin black rings. We stress, however, that our construction is not the most general 
one to 2nd order in the derivative expansion. We will be ignoring terms that are constructed from the 
Levi-Civita symbol and which are dimension-dependent. The reason for this will be explained in Sec. 5. 

This paper is organized as follows. In Sec. 2 we lay down the general framework for the variational 
calculus and the construction of the free energy in a derivative expansion. We then show how to iteratively 
account for higher order derivative corrections and obtain expressions that exhibit the coupling between 



^Note that in the case of spinning corrections to the motion of point particles, it was aheady known that the equations of 
motion were captured by a multipole expansion of the stress-energy tensor, see for example [44]. In the context of curvature 
corrections to cosmic strings, this was also derived in [45]. Here, we have generalized these results to arbitrary space-time 
dimension and for p-branes. 



the hydrodynamic and elastic modes. Sec. 3 is quite technical but equally important. There we show how 
these corrections can be understood as a multipole expansion of the stress-energy tensor and proceed to 
construct conserved currents and charges. In Sec. 4 we match our results with the measurement of the 
Young modulus for black branes, firmly establishing what is meant by an elastic expansion for black holes. 
In Sec. 5 we comment on open issues and future work. We also include App. A where we describe in 
detail the notation used in this paper and collect different results of variational calculus and of geometry 
of embeddings. In App. B we analyze the boundary conditions for hydrodynamic modes. 

2 Effective action for (non)-extrenial branes 

In this section we construct the most general free energy for a neutral stationary fluid brane quadratic in 
the extrinsic curvature and extrinsic twist potential to second order in a derivative expansion. We consider 
a (p + l)-dimensional surface with worldvolume topology Wp+i = M x i3p embedded in a background 
space-time of -D = n+p + 3 dimensions with metric Qfiuix'^) and space-time coordinates x". The position 
of the surface is described by a set of mapping functions X^{a°'), where the set of coordinates a"" cover the 
embedded submanifold. The indices fi, u, A... are space-time indices, the indices a, b, c... denote directions 
along the worldvolume surface while the indices i,j, k... denote transverse directions to the worldvolume. 
The worldvolume inherits an induced metric of the form ^ab = g^iuUaWj^, where Ua = daX^, and if 
the worldvolume is bent one can assign to it an extrinsic curvature defined as Kab^ = ^ aU^ , where the 
covariant derivative Va is defined in App. A. Any space-time vector u^ can be decomposed into tangential 
and orthogonal components to the worldvolume surface such that v^ = v"'Ua + v^n^, where the set of 
tangential and orthogonal projectors satisfy gf^u^au'^ = 0. Given the set of orthogonal projectors one can 
define the extrinsic twist potential uia^ = —n^^Van^^, which is anti-symmetric its two transverse indices, 
as well as the outer curvature O.ab'^ associated with it (see App. A). It is also useful to introduce projectors 
onto the worldvolume and onto its transverse space while keeping track of space-time indices. For that 
matter we introduce the first fundamental tensor 7^^"^ defined as 'y^'^ = ^"'^u'^au'^b , which projects onto 
Wp+i , and the orthogonal projector _L'^'^ given by the expression J-^'^ = g^^ — ^^^ . Further, we assume 
the existence of a set of commuting worldvolume Killing vector fields that we take to be of the generic 
form^ : 



k»9a = dr + Yl ^^"^^'^'^ ' ^" = (^' <^- ••■) ■ (2-1) 



a=l 



Here the set of angular velocities OS'^' is constant. We begin by reviewing the 0th order action and 
generalizing the analysis of [21] while simultaneously highlighting some of the elastic properties of fluid 
branes. The method for iteratively correcting the action then follows. 



■^The existence of such worldvolume Killing vector fields is a requirement for stationarity of the overall configuration (see 
for example [52]). 



2.1 0th order fluid-elastic system 

To 0th order in a derivative expansion the surface is described by an induced metric ■jab and the set of 
commuting Kilhng vector fields (2.1). The only natural scalar invariant that can be constructed from 
these two is the modulus of the Killing vector defined as k = | — 7a(,k°k p^. Therefore, the action must 
be a functional of k^, 

T[X^^] = f C {V^, k) = / V^ Ao(k) . (2.2) 

JWp+i jyVp+i 

Our goal is to make a small deformation 5X^ of the worldvolume geometry along both tangential and 

orthogonal directions. To this aim we decompose the deformation as 6X^^ = ^'^ = ^°'Ua + $*n^. Under 

arbitrary small variations the metric changes by a Lie derivative such that 

Slab = 2V(,$,) - 2Kab'^i . (2.3) 

Further, we assume the set of angular velocities in (2.1) to be held constant during the variation and 
hence the variations of the Killing vector fields 6k^ are zero. Therefore, the variation of the action (2.2) 
can be written only in terms of the variation of the induced metric. This fact leads to a result of the 
form: 



6T[X>'] = / V^ - V J T'^^'^b + ^bVaT'"' + T'^'Kab''^, , (2.4) 

JWp+i \ V / / 

where we have defined the monopole source of stress-energy tensor in the usual way, 

T-' = -^^ • (2.5) 

Before we analyze the explicit form of (2.5) note that in order to obtain the equations of motion from 
(2.4) one must require 5J^[X^] = 0. The first term in (2.4) is a boundary term and can be written in the 
form 

- [ ^/^T'^Sa^b , (2.6) 

JdWp+i 

where r]"" is a unit normal vector to the brane boundary and h the determinant of the induced metric on 

the boundary. Hence, a well posed variational principle requires T°'^rja\dv\! +1 = 0. The second term in 

(2.4) results in an identity which is trivially satisfied as we will show below, 

VaT"^ = . (2.7) 

In total, the variational principle (2.4) yields the non-trivial equation of motion 

T^^Kab' = . (2.8) 

Eq. (2.7) expresses stress-energy tensor conservation and is responsible for the fluid behavior along world- 
volume directions while Eq. (2.8) contains the elastic degrees of freedom which are manifested along trans- 
verse directions to the worldvolume. These equations will be analyzed in greater detail as we progress in 
this section. 



For explicit examples of non-extremal brane actions see [19, 53, 54, 55, 56]. 



The explicit form of the stress-energy tensor 

To evaluate the stress-energy tensor (2.5) we need the variation of k with respect to 7^6, this is simply 

5k = -|nV<57a6, ^" = Y- (^-^^ 

For an action of the type (2.2) the stress-energy tensor takes the following form: 

j.ab ^ yab ^ _Ao(k)7'^^ + A^(k)knV . (2.10) 

Here Ao(k) indicates a derivative with respect to k. This is easily recognized as the stress-energy tensor of 
a perfect fluid and indeed leads us to identify the set of normalized vectors u^-Ua = — 1 as fluid velocities, 
while thermodynamics allows us to identify the pressure P and energy density e such that 

P = -Ao(k) , e + P = A'o(k)k. (2.11) 

Assuming the Gibbs-Duhem relation e + P = Ts and knowing that for such systems the constant global 
temperature T is related to the local temperature via T = kT [52], using the identification (2.11) one 
obtains an expression for the local entropy density 

. = ^A;,(k)k2 . (2.12) 

Given these identifications, the set of Eqs. (2.7) are interpreted as the fluid equations on the worldvolume 
Wp+i. However, the stress-energy tensor (2.10) satisfles the conservation equation (2.7) regardless of any 
thermodynamic interpretation. Indeed, using (2.10) in Eq. (2.7) results in 

-V%(k) + J^A;,(k)V^k2 = 0, (2.13) 

which is trivially satisfied. In order to obtain the above, one needs to use the fact that the expansion 
6 = Vflii" vanishes by virtue of the Killing equation Vak;, = Vr^ky and that for any worldvolume scalar 
or tensor T, the Lie derivative along the Killing vector field vanishes X^T = 0. 

The elasticity of fluid (nieni)-branes 

Eq. (2.8) encodes the elastic degrees of freedom of the brane. In fact, note that it has a direct analog in 
terms of classical elasticity theory: it is the equation of motion obtained by varying the free energy of a 
thin stretched membrane when bending effects can be neglected [57, 21]. For very thin membranes, the 
effect of bending is always subleading when compared to deformations caused by stretching or compression 
[57]. To further see the connection with elastodynamics it is instructive to imagine the following scenario. 
Suppose that to a given fluid conflguration satisfying Eq. (2.8) one applies a small deformation of the 
embedding surface in an arbitrary orthogonal direction <&*. Prior to the deformation the metric jab 
measured distances between fluid elements on the surface, while after the deformation the actual value 



of 7a(, measures the new distances on Wp+i. Therefore 7ab describes the state of strain of the brane and 
one can define the Lagrangian strain tensor Uab"^ as [21] 

Uab = -^{lab-%b) . (2.14) 

For infinitesimal displacements along <!>*, the strain tensor changes by a Lie derivative and hence 

£<i,'Uab = dUab = --d-fab = AV^, . (2.15) 

Thus, we conclude that the extrinsic curvature tensor Kab^ measures the strain induced on the brane 
due to a deformation of the surface along orthogonal directions [5, 21]. Making the extra assumption 
of the existence of a background Killing vector field k'* whose pullback onto the worldvolume coincides 
with the worldvolume Killing vector field k" allows us to write Tsn^ = —sdfjj'. This, together with the 
form of the stress-energy tensor (2.10) and the identification (2.11), allows us to rewrite Eq. (2.8) when 
contracted with the deformation vector $* as 

dP = -PdV , (2.16) 

where we have defined the relative change in volume along an orthogonal direction as dV = {l/2)^°'^d'yab- 
Eq. (2.16) allows for the definition of the modulus of hydrostatic compression /C that measures the brane 
response to variations in volume such that 

1 fdV\ 1 



(2 \7\ 
K ydPj^ Ao(k) ■ ^ ' ' 

Further elastic properties can be highlighted by using the first law of thermodynamics for the fiuid 
de = Tds and defining the solid density p = e + P = AQ(k)k which along the directions $* expresses the 
first law of thermodynamics for an elastic membrane that has been subject to hydrostatic compression 
dp = Tds — PdV [21]. To every fluid membrane one can also assign an elasticity tensor that measures the 
deformation of the stress-energy tensor (2.10) due to surface deformations. In order to see this precisely 
one performs a small deformation of the equations of motion (2.8) along the vector $*. To this aim, one 
requires the deformation of the extrinsic curvature tensor along orthogonal directions (see App. A), 

5^Kj = n\Va (±'^.V6($%^,)) - R'baj^' - Kac'K%^^ , (2.18) 

where R^uXp is the Riemann curvature tensor of the background geometry. Using this transformation 
rule, an infinitesimal deformation of Eq. (2.8) yields 

E-b'^dj^jKj'^j + r«''nV.V„V,(<l>^n'^,) = T''''R\bj^' , (2.19) 

where we have defined the elasticity tensor for fluid branes E°'^'^'^ through the Hookean relation dT""^ = 
E°'^'^'^dUcd, which for systems of the type (2.2) takes the general form 

E'^^'^'' = -2 f Ao(k)7'^('=7')' - f ^^1 1^" - 2 f ^^l") . (2.20) 



This is the usual definition of Lagrangian strain [58] but now including the timelike direction as well. 



The structure of (2.20) is that of a material characterized by a varying moduhis of compression. The 
case of Dirac-branes is obtained when Ao(k) is constant and the two last terms vanish. The fact that 
the elasticity tensor is only probing compression and stretching is because we are working in the limit in 
which the fluid is confined to an infinitely thin surface as it will become clear in Sec. 3. We note that 
the linearized equation (2.19) has been previously obtained by Carter in [37], here we merely applied it 
to a fluid brane. It is interesting to note that Eq. (2.19) already encodes some of the structure of the 
equations of motion when bending effects are taken into account. 

The free energy interpretation 

The 0th order action (2.2) can be interpreted as the free energy of the fluid-elastic system [19, 54]. After 
identifying the stress-energy tensor (2.10) with that of a perfect fluid, then by defining the local Gibbs 
free energy density Q and using (2.11) one obtains 

g = e-Ts = Xo{k) . (2.21) 

The integral over the worldvolume of the density Q yields the total free energy of the system (2.2). This 
interpretation can also be realized at the level of the global charges. First note that the first term in 
the variational principle (2.4) is a total divergence term. For the action (2.2) to be invariant under an 
infinitesimal shift of worldvolume coordinates along a worldvolume Killing vector field, one must have 
that 

Va (r^^i^b) = . (2.22) 

This is satisfied due to the Killing equation and the symmetry of T""^. The set of surface currents T"*k{, is 
conserved and with those it is possible to construct a set of conserved surface charges (mass and angular 
momenta) of the form 

M= [ dYr^.T-'naCb , J^"^ = - I dV,p.T''\axt^ ■ (2.23) 



\p)-L nal;b , J'-'=- aV^p)l UaXb 

J lip JBp 

In writing these expressions we have assumed that the worldvolume timelike Killing vector field ^""da = d-r, 
whose norm is the redshift factor Rq, is hypersurface orthogonal with respect to Wp+i. Further, we have 
introduced the spatial measure dVip) on the worldvolume, a unit normal vector n"" = S,'^/Ro orthogonal to 
spacelike slices of Wp+i and defined the rotational Killing vector fields in (2.1) as x da = dfpa. In Sec. 3 
we will show that this set of surface charges is indeed the set of global conserved charges associated with 
the fluid brane. Using the definition of the local entropy density (2.12) it is straightforward to construct 
the global entropy of the system from the entropy current J° = su"": 



S = - dV^p^J^a ■ (2.24) 

JBp 

Having defined the charges (2.23) and the global entropy it is possible to verify that T = M—Q^°-\j(^a) — TS 
and hence the variational principle dJ^ = with fixed n'"-* and global temperature T, requires the first 



law of thermodynamics to be satisfied [19, 54]: 

dM = n^^Ujf^a) + TdS . (2.25) 

2.2 General framework for higher order corrections 

We will now show how some of the ideas of the previous section can be pushed to second order in a 
derivative expansion. The general method, following [46, 50] , consists in adding all possible scalars 
constructed out of derivatives of worldvolume quantities to the action (2.2). We will split the type of 
corrections that can be supplemented to the action in three parts. Hydrodynamic corrections are those 
which only involve derivatives of fields that characterize the intrinsic geometry of the brane. These are, 
for example, terms proportional to Vak, Vjaku, to the worldvolume Riemann tensor TZabcd or to the 
tangential projection of the background Riemann tensor Rabcd- Elastic corrections, on the other hand, 
are scalars proportional to the extrinsic curvature tensor Kab^. Note that by (2.3) or (2.14), the extrinsic 
curvature tensor is a one-derivative term along transverse directions to the worldvolume. Finally, spin 
corrections are corrections proportional to the extrinsic twist potential uja^ , which is also a one-derivative 
term. This means that we will consider a generic free energy of the form 

T[X>']= I £(V=7,7a6,k«,V„K,bSL^,*^). (2.26) 

Since the Killing vectors k" are held constant during the variation, variations of (2.26) can be analyzed 
solely through variations of the metric 7af,, the extrinsic curvature tensor Kab and the extrinsic twist 
potential oJaK For this purpose, we define the dipole moment P « and the spin current S°'ij of the 
fluid-elastic system as^ 

V-\ = ^^^ , 5% = ^^ . (2.27) 

To study small deformations of the geometry (2.26) along both tangential and orthogonal directions 
to the worldvolume it is required the knowledge of the variation of Kab as well as of ooa^ along these 
directions. Using (2.18) for the transverse variations of Kab , we have that for the tangential variations 

5\\Kab' = ^VpV.Kab" + 2i^,(,^V,)$'= , (2.28) 

while for the extrinsic twist potential we have that (see App. A) 

6u;J' = nVn^^VfeC^/^ + cvb^'Va'^'' - 2Kab^'n^^pV\'^''nkP) + R'ha1>^ , (2.29) 



^The reason for the interpretation of the quantities (2.27) as the dipole moment and spin current wiU become clear in 
Sec. 3. 



where w^'^^ = n^inP jUOaK Using the transformation rules (2. 18), (2. 28) and (2.29), the total variation of 
(2.26) takes the form 

Jwp+i L ^ ^ 

where we have used the definition of the outer curvature ^ab^'' (see App. A). As in the 0th order case 
(2.4), the variation yields a total divergence which can be integrated to a boundary term. However 
note that any component of the form Vfe<&^ is not independent on the brane boundary, indeed one can 
decompose it such that Vh<l>'^ = r]i,7]"''V a^^ + vi,"''Va^P, where Va'^ are boundary coordinate vectors and 
the indices a,b, ... label boundary directions. The normal component to the brane boundary i]'^'S/a^^ is 
independent but Va^Va^P is not. Hence, assuming that the variations ^^ vanish on the boundary of the 
brane boundary itself (if existing), a well defined variational principle requires 



V'^'^VamlaWp+r = , cS'^*^??a|0W,+ i = 



Va (v^'^n^iTjaVt^) + rja (r^'ub^ - V^'^Kj^'ub^ + ^^V^P'^^'' - 2n^jS'^'Kb^, 



(2.31) 



law, 



p+i 



. 



Further, the second term in (2.30) yields a set of non-trivial identities that must vanish. With the help 
of the Codazzi-Mainardi equation 

R'cba = n'p{VbKa/-VaK,bn , (2.32) 

and the Ricci integrability condition [59, 38] 

Rab''=nJ^-KjKb'^ + KbjKa'^ , (2.33) 

it can be brought to the form 

VaT'^' + u^'VaVcV'''^ - 2 S^jKajK'^^ = V^'''R\,i - S'''m\ij + J'^VaS\j . (2.34) 

In Sees. 2.4 and 2.5 we will give a few examples of how this equation is satisfied for the actions we consider. 
Eq. (2.34) can be seen as the modified intrinsic dynamics of Eq. (2.7). The reader may wonder if there 
is another definition of the stress-energy tensor for which Eq. (2.34) would just express its conservation. 
Indeed in Sec. 3 we will show that it is possible to define the linear momentum which is always conserved 
in flat space. Finally, the variational principle (5J-"[X'^] = yields the non-trivial equation of motion: 

T'^^'Kab' = -n'pVaVbV^'P + 2 n'pVb (SaP'K'^'j) - V^'' R\jb - S'''^^ R\kj ■ (2.35) 



This is of course the modified version of the equation of motion (2.8) due to the presence of dipole 
jyabi g^j^^ spinning S°'ij effects and, in fact, it is the generalization to arbitrary co-dimension and to 
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curved backgrounds of the classical equation of motion of a deformed thin membrane when bending 
effects and rotation in transverse directions are taken into account. To see this precisely, let us focus on 
a flat background and consider the well studied case of the Helfrich-Canham bending energy (1.1) for 
co-dimension- 1 surfaces and in the absence of spinning effects, that is, S'^ij = 0. In this case, as we will 
show below, we have P'*''* = 207"* i^T* where K^ = ^""^K^b is the mean extrinsic curvature vector. The 
equation of motion (2.35) can then be rewritten as 

2 a n\V^XP + V^Kab' = , (2.36) 

with V = T — 2aK iK^ and where V^ is the square of the Laplacian, or the biharmonic operator 
usually found in classical elasticity theory [57]^. Therefore, the theory described by Eqs. (2.31)-(2.35) is 
a general relativistic generalization of classical elastodynamics of thin membranes. 

It is also possible to construct surface currents and charges associated with the systems (2.26), how- 
ever, since there are subtleties that can only be resolved when relating this approach to a multipole 
expansion of the stress-energy tensor, we will leave that for Sec. 3. On the other hand, the construction 
of the global entropy and entropy current as in (2.24) to second order in the derivative expansion will be 
lacking in this work as it has not been developed for the systems we consider. Such endeavor is possible 
to accomplish and it will be presented in a future publication [60] . 

For now, we will study the different contributions that can arise at each order in the derivative 
expansion. To each term involving one derivative we associate it with the expansion parameter e. Two- 
derivative terms are of order O (e^)- For clarity of explanation, we decompose the possible corrections 
to the intrinsic stress-energy tensor T"'^ into hydrodynamic 11"^, elastic r"'' and spin 0"'' contributions 
such that 

rj^ab ^ j.afe ^ jjafe ^ ^ab ^ Qab ^ ^2.37) 

The higher order corrections 11'^*, r"^, 0"^ as well as the dipole moment 2?°^* and the spin current S°'ij 
will be the sum of the contributions from the different scalars that can be added to the action (2.26). 
Generically, 

n"^^ = "^uf , t"^ = ^t^\ q'^^ = ^ e^^ v^ = Y^ v^^' , s""'^ = Y^ s^'^ . (2.38) 

a a a a a 

Note however that the hydrodynamic corrections do not contribute to the dipole moment T) neither to 
the spin current S^^^ by virtue of the definition (2.27). 



®Note that in classical elasticity theory V^'' is conserved while here, due to Eq. (2.34) it is not. However, the tensor 
pan ^ po^u^M _|_ ^^VaC" is conservcd in flat space and V'^Kap,' = V^Kab'- Note also that for co-dimension- 1 surfaces, 
the extrinsic twist oja'"' vanishes. 
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2.3 1st order action 

To 1st order in the derivative expansion the only possible terms that can be added are of hydrodynamic 
nature. This would be terms of the form: 

k'^Vak , Vak'^ . (2.39) 

However these terms vanish because k*^ is a Killing vector field. This is in agreement with the analysis of 
[46, 50] for both stationary and non-dissipative fluids. Indeed, for a generic dissipative fluid the entropy 
current to first order in derivatives is proportional to the square of the expansion 9 = Vau"" and to the 
square of the shear a . These two quantities vanish for stationary flows. One may also consider extrinsic 
curvature corrections to 1st order in derivatives but due to the transverse index in Kab^, any scalar built of 
the extrinsic curvature needs to be contracted with another copy of itself. Therefore, extrinsic curvature 
corrections only enter to second order in the derivative expansion. The same argument holds for terms 
proportional to the extrinsic twist potential^. This is in agreement with the results obtained from gravity 
for the leading order corrections to thin black rings [17] and blackfolds [16] for co-dimension surfaces 
higher than three for which a system of the type (2.26) is supposed to be the correct description [5]. An 
exception to this argument of absence of extrinsic corrections to first order is the case of co-dimension- 1 
surfaces where we can simply omit the transverse index from Kat^. These cases will be dealt with in 
Sec. 2.7. All in all, for surfaces of co-dimension greater than one we conclude that there are no corrections 
to first order in the derivative expansion. 

2.4 2nd order elastic corrections 

To second order in the derivative expansion, focusing on the elastic corrections, we can add a total of 5 
different terms to the action (2.26) which are quadratic in the extrinsic curvature. These terms are of 
the following form: 

(2.40) 

A4(k)k'^k^i^,6^i^, , X5{k)k''k\'^k''Kab'K,a^ . 

The first term in (2.40) is the generalization of the Helfrich-Canham bending energy (1.1) for surfaces of 
arbitrary co-dimension and for a non-trivial response coefficient Ai(k). The first two terms in (2.40) for 
the case in which both Ai(k) and A2(k) are constant build up Carter's model of stiff strings and branes 
[36] and have also been extensively studied by Guven [38, 61]. The other three terms have not been 



^The shear tensor is defined as a"'' — p'^'^p'"^ I V(citd) ~ -"fed ) , where P"' is the projector along transverse directions to 

the fluid flows and defined as P"* = 7"'' + u^'u''. 

*If we considered terms constructed from the Levi-Civita symbol then it would be possible to add a term of the form 
k^EijUa'''- However, even though these are not considered in this paper, they are of physical interest, as they may be useful 
in the description of spinning black holes (see Sec. 5). 
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considered previously in the literature and they constitute the generalization of the Helfrich-Canham 
bending energy for stationary fluid (mem)-branes. We note that the reason why the second term in 
(2.40) was not necessary to be added to the action of cellular membranes is that the first two terms are 
coupled to each other and in certain conditions, such as in flat space, they can be shown to be equivalent. 
We will explain this in detail in Sec. 2.6 when we study the coupling between modes. Decomposing each 
scalar in (2.40) as XaO^)jC-a, we summarize below the contributions of each term to t and 'D^'"^: 



Scalar 


a 


'Tjabi 
a 


Ai(k)£i 


-Ai(k)£i7"^ + A'i(k)k£i'uV + AXi{k)K''\K' 


2Ai(k)7'^*if^ 


A2(k)£2 


-A2(k)£27"'' + A'2(k)k£2w"?^^ + 4A2(k)A'''^iK^* 


2A2(k)i^'^''* 


H^)C3 


-A3(k)£37"'' + X^{k)kCsu''u' + 2Xs{k)k^k'' K\iK''d 


2A3(k)k'^k(»i^'')rf* 


A4(k)£4 


-A4(k)£47'^^ + X'^{k)kC^u''u'' + 2A4(k)k^k'^i^'^''ii^cd^ 


A4(k)k'^k^K' + A4(k)7'^^k^k'^i^cd^ 


A5(k)£5 


-A5(k)£57''' + A'5(k)k£5n'^n^ 


2A5(k)k"k^k^k'^i^,d^ 



We see that t""^ is the sum of a total of 5 contributions t"''' = "^^=1 '^a' ^^"^ contains a perfect fluid part 
(see (2.10)) and an elastic part. When r"^ is added to the 0th order contribution (2.10) we decompose it 
into the total perfect fluid part T"^ and elastic part such that 



where the elastic deformation to the worldvolume stress-energy tensor S"-"'^"- has the form 
gabcd ^ 2 ('2Ai(k)7'^S"^ + 2A2(k)7'^'^7''^ + A3(k)7'^'^k^k'= + A4(k)7"''k'=k'^^ 



(2.41) 



(2.42) 

Rewriting the 0th order contribution (2.2) as Ao(k)£o with Cq = I and using the identification (2.11), 
the total fluid part T"^ can be written as 



o-„,b 



where 



7-ab ^ p^ab + (^ + V)u''lJi 



V = -Y,^c.{\^)^» , £ + V = Y,X'^{k)kCo 



(2.43) 



(2.44) 



«=o 



«=0 



Indeed, the quantity defined above as the modified pressure V is nothing more than the sum of all scalar 
contributions multiplied by a minus sign, that is, the 0th order scalar (2.2) plus the 2nd order ones (2.40). 



^Note that elastic corrections do not contribute to the spin current S"-' . 
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The dipole moment 'D°'^^ is also the sum of 5 different contributions P'^''* = "^^^^ T)^^ and can be put 
into the elegant form 

pafci ^ y^bcdj^j ^ ^2.45) 

where we have defined the Young modulus of the fluid brane as 

yabcd ^ 2 (x^j^^j"'^ + AaT*"^ V^^ + A3k("7'')(^k'^) + — (j''\^k'^ + 7'="'k''k^) + Agk^'k^k^k"'') . (2.46) 

Here we have omitted the dependence of the response coefficients Aq on k. The Young modulus (2.46) 
is the physical quantity that encodes all the possible responses of stationary fluid branes to bending 
deformations. Furthermore, it exhibits all the symmetries of the usual elasticity tensor of a classical 
anisotropic crystal y"-^'^'^ = y{°-b)i<^d) _ ycdab ^ ^g ^^ ^jj^ g^g -^^ g^^^ ^^ ^j^g Young modulus measured from 

gravity for bent black branes is a particular case of (2.46). Given the structures (2.45)-(2.46) and the 
identification (2.11) we present two useful and equivalent rewritings of the action (2.26) including only 
elastic corrections: 

n^^] = - I V^ f P(k) - Iv^^KaA =- f V^ f P(k) - hab^'K'^^K.A . (2.47) 
JWp+i V ^ / JWp+i \ ^ J 

In Sec. 4, we will study an example of this action. However we note that this action is exactly the type 

of action expected for the bending of thin membranes from classical elasticity theory [57], here we have 

presented a relativistic generalization of it. 

Vanishing of the intrinsic equation of motion 

In Sec. 2.2 we mentioned that the intrinsic equation of motion (2.34) resulted in mere non-trivial identities 
for the actions we consider. One can verify that this is indeed the case for all corrections presented in 
the table above. Here we will just present two examples of how this is done. For the mode Ai(k), when 
introducing rf'' and T>'f''^ into Eq. (2.34) a part of the stress-energy tensor is conserved due to (2.13), the 
remaining part is brought to the form 

2Ai (k)i^i iv^'K' - VaK^'^') = 2Ai (k)7"=K*i?^,e • (2.48) 

Using the Codazzi-Mainardi equation (2.32), the l.h.s. of the equation above can be seen to be equivalent 
to the r.h.s.. In fact, for the terms Ai(k) and A2(k) one only needs to use (2.13) and (2.32) to show the 
identity (2.34). We now take the case of the term A5(k). After using Eqs.(2.13) and (2.32) we are left 
with 

As (k)k'^k^k^i^,/i^ae^ V^k'^ + As (k)k» Va (k^k'^k'^Ke/i^c^) = . (2.49) 



Now we need to remember that the Lie derivative along any worldvolume Killing vector field of any 
world volume tensor must vanish, that is, 

£^ (k^k^k^Ke/i^c^) = k-^V, (k^k'^k'^Ke//^c'i) - k^k'^k^^Ke/Z^ac.V'^k'' = . (2.50) 

When Eq. (2.50) is used, Eq. (2.49) is trivially satisfied. This fact, together with the Killing equation, is 
enough to verify the identity (2.34) for the terms A3(k) and A4(k). 
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2.5 2nd order spin corrections 

Spin corrections to O (e^) in the derivative expansion are the terms that can be added to the action 
which are quadratic in the extrinsic twist potential. These do not contribute to the dipole moment 2?"''*. 
Therefore, ignoring terms constructed from the Levi-Civita symbol, we can add two different scalars to 
the action 



Wi(k)uj"-ijUJa^ , tZ72(k)k''k OJaijUJb 



«J 



(2.51) 



Denoting each term by Wa{'k)Wa we summarize below their contribution to the stress energy tensor G'^'' 
and to the spin current S*^^^: 



Scalar 


ef 


S^'' 


t^i(k)>Vi 


-roi(k)Wi7"^ + ro;(k)kWin«n'' + 2wi{'k)uj'''iuj'',j 


2roi(k)a;"*J 


UJ2ik)W2 


-ro2(k)W27''^ + CT^(k)kW2u"u'' 


2w2{k)k''k^ujb'^ 



The stress-energy tensor of these contributions contains a perfect fluid part that can be incorporated in 
the total pressure defined in (2.44). Defining a spin deformation tensor of the form 



^ab 



ai,6 



2roi(k)7"'' + 2tz72(k)k''k'' , 



(2.52) 



the action (2.26) including only the 0th order contribution (2.2) and spin corrections can be written in 
two useful ways: 



-^[^1 = - / V^ f ^(k) - ls'^^JU^aA =- [ V^ (P{k) - ^S'^'ua'^Ub.j) . (2. 

We will see however that the term wi (k) is not a desirable one in a physical theory. 



53) 



Vanishing of the conservation equations 

As in Sec. 2.4 the vanishing of the conservation equation (2.34) can be shown by following the same 
steps as in the previous section together with the definition of the outer curvature (see App. A) and the 
Ricci integrability condition (2.33). However, there is a subtlety in this demonstration. For the scalar 
ti72(k)W2, one can show that the last term in Eq. (2.34) vanishes independently, that is. 



n'n^xVaS^"^ = 



(2.54) 



a result that follows from the fact that £k (k w;,*-') = \i.""n^ pU^ \S/ a (k Wf,'' 



c i,!iP^\ = 0. However, even though 
the scalar ■c<7i(k)Wi satisfies (2.34), it does not satisfy (2.54) individually. Eq. (2.54) expresses conserva- 
tion of angular momenta in transverse directions to the worldvolume, as we will see in Sec. 3. Therefore 
the term ■ci7i(k)Wi is explicitly violating this conservation. 
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2.6 2nd order hydrodynamic corrections and mode coupling 

In this section we examine the possible 2nd order hydrodynamic corrections to the action (2.26) and the 
stress-energy tensor T""^ as well as the coupling between these modes and the ones found in the previous 
sections. For all these corrections the dipole moment P vanishes. There are 7 types of scalars that can 
be constructed from the intrinsic geometry^'': 

t;i(k)VaV'^k , U2(k)7^ , V3ik)k^k'TZab , 

(2.55) 

U4(k)V[,k,]V["k''] , t;5(k)V,kV'^k , Veik)R\a' , V7ik)k''k'R'^acb ■ 



Here we have introduced the Ricci tensor TZab and Ricci scalar TZ of the worldvolume geometry. The 
terms f6(k) and f7(k) can be related to some of the other hydrodynamic terms and to the elastic modes 
in (2.40) through the Gauss-Codazzi equation as we will explain below. This of course expresses the extra 
difficulty of the problem compared to [46, 50] since not only couplings to the worldvolume geometry need 
to be considered, but also couplings to the background geometry. Further, the terms t'i(k), f3(k), f4(k) 
and v^Ck) are not independent of each other and are related by on-shell equivalences. To see this note 
that 



/ V^V3{k)k^k'nab = - [ V^Vsik) [Va, Vb] k"k^ 

= I V^V3{k) (VakV^k + kVaV^k - V[,kb]V["k^] 



(2.56) 



Moreover, the modes ui(k) and v^ik) are directly related to each other, modulo a boundary term, such 
that, 



/ 



7t;i(k)VaV"k = - / V^^^i(k)VakV''k + boundary term . (2.57) 



p+i -J yvp+i 



Therefore, using Eqs. (2.56)-(2.57) we can eliminate the terms f4(k) and v^{k) and only consider the 
first three in (2.55) ^^. The scalars (2.55) yield boundary contributions that were not included in the 
general analysis of (2.30). For this reason we have prepared App. B where these issues are clarified. 

In order to perform orthogonal variations of (2.55) it is only necessary to know how the worldvolume 
Christoffel symbols ^ab'^ vary with respect to the induced metric (see App. A). Denoting each term by 
Va{k)Va^ we summarize below the contribution of each of these 3 relevant hydrodynamic terms to the 



^"Remember that since we are dealing with stationary flows both the expansion 6 and shear a"^^ are zero and hence any 
term proportional to these vanishes. These are, however, non-zero for generic non-dissipative flows [50]. Further, the reader 
may wonder why we have not considered couplings of the form R^-' ij or contractions with Waib- The reason is that these 
terms are not well defined on the worldvolume and are related to geometric quantities with support on the transverse space. 
See [62] for a setting in which these terms play a role. 

^^This is exactly the same type of analysis as in [50] but now applied to a submanifold embedded in a background 
space-time. 
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stress-energy tensor IT : 



Scalar 




viO^Wi 


-r;i(k)Vi7'''' + v[{k)kViu''u'' + -/"''Vc (i;i(k)V^k) + kn'^n^VcV^t;i(k) + 2t;i(k)V'^V^k - 

2V('^(t;i(k)V^)k) 


t;2(k)V2 


-U2(k)V27"'' + i;^(k)kV2?/"n'' + 2^;2(k)7^«'' - 2V'^V''i;2(k) + 27«^VcV^i;2(k) 


t^3(k)V3 


-U3(k)V37''' + v'^ik)kV3u''u'' + VeV^ (t;3(k)k'^k^) - 7'^'VeVrf (i;3(k)k^k'^) + 

2VM^ (v3{k)k''^k''') 



Looking at the table above, we see that again there are contributions of the perfect fluid form which can 
be included in the definition of V and £ in (2.44), and other new contributions which are second order 
in derivatives of the intrinsic variables. It is also possible to show that all such contributions satisfy the 
identity (2.34). The only necessary ingredients to show this are Eq. (2.13), the Codazzi-Mainardi equation 
(2.32), the fact that the worldvolume Einstein tensor is divergenceless VaG''^ = Va [Tl"^ - (l/2)7''*7^) = 
and stationarity of the overall configuration. 



Coupling between elastic and hydrodynamic modes 

Here we study the coupling between the elastic modes (2.40) and the hydrodynamic modes (2.55). The 
most important piece of knowledge is the Gauss-Codazzi equation [59], 



Rabcd = T^abcd — KacK^i + Kad^Kbc 



(2.58) 



relating the intrinsic geometry to the extrinsic geometry. We begin by contracting (2.58) with the induced 
metric 7^"^, this leads to 



7 Rabcd = Tlac - KacKi + KabK c 



(2.59) 



If further contracted with k^'k'^, this equation expresses the coupling between the hydrodynamic terms 
V2,{k),VT{k) and the two elastic contributions A3(k),A4(k). Indeed we can see why we can eliminate 
V'jik) in favor of the other three terms. Hence, if ^^'^Rabcd is non-zero, only 3 of these 4 terms can be 
regarded as independent, while if TZac or ^ Rabcd is vanishing, only 2 of the 3 remaining terms can be 
regarded as independent. However, if both the background and the worldvolume are flat, only 1 of the 
elastic contributions is independent. In this case, the two elastic terms left are exactly equivalent to each 
other. If we now contract (2.59) with 7"^^ we obtain 



,,ac^hd 



ab 



r^Y'^Rabcd = n- K'Ki + K^'^Kab 



(2.60) 
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This equation, in turn, expresses the coupHng between the hydrodynamic terms vi{]i.) , VQ{\i) and the 
two elastic constributions Ai(k), A2(k). Therefore, if the background is curved, only 3 of these 4 terms 
can be regarded as independent and again, if R"'ba or 7^ vanish, only 2 of the 3 remaining terms are 
independent. If both the background and the worldvolume are flat only 1 of the 2 elastic contributions 
is independent. For the fluid membranes studied in theoretical biology for which the Helfrich-Canham 
bending energy (1.1) was crucial to understand, the term f2(k) = constant, even though entering at the 
same order in derivatives, is not added to the action as it is a purely topological invariant. Moreover, in 
flat space, and since the two elastic contributions Ai(k) and A2(k) are coupled through the relation (2.58), 
there was no need to consider A2(k) in describing red blood cells. On the other hand, the hydrodynamic 
contribution f i(k) and the elastic contribution A5(k) are independent and can be excited without exciting 
the remaining ones. Furthermore, the spin modes (2.51) are decoupled from the hydrodynamic and elastic 
modes to order O (e^)- 

Summarizing our results, for surfaces of co-dimension greater than one, if the background and the 
worldvolume are curved, the fluid-elastic system (2.26) admits a 9 parameter family of stationary fluid 
branes^^. On the other hand, if either the background or the worldvolume are flat, there exists a 7 
parameter family of stationary fluid branes. Finally, if both the background and the worldvolume are flat, 
there exists a 5 parameter family. The presence of elastic degrees freedom introduces further non-trivial 
response coefficients than those analyzed in [46, 50] for both stationary and non-dissipative space-filling 
fluid flows. 

2.7 Fluid membranes: co-dimension- 1 surfaces 

As mentioned in Sec. 2.3, in the case of co-dimension- 1 surfaces there are non-trivial elastic corrections 
to 1st order in the derivative expansion. Since the transverse space is only 1-dimensional, the extrinsic 
curvature tensor can be written as X^ = Kab ■ We omit the transverse index and only use it if necessary. 
For this reason, there are two non-trivial and independent elastic contributions to 1st order in derivatives 

^i(k)i^ , ^2(k)k'^k''i^,b . (2.61) 

These terms contribute to T"^ and V"-^-^ while the equations of motion (2.31)-(2.35) remain the same but 
have only one component along the transverse space. The first term in (2.61) had been previously con- 
sidered in the literature [38] but it does not play a role in real fluid membranes since the mean extrinsic 
curvature is required to be kept constant due to experimental constraints, that is, the variation of (2.26) 
must be supplemented with local constraints [39]. The second term in (2.61) has not been previously 
considered in the literature and is accounting for the fact that the fluid may be rotating. Denoting each 
scalar in (2.61) by '&aO^)Ca and their contributions to T"^ as f"^, we summarize the results below: 



^We are taking into account the contribution zu2(k) but not ci7i(k) since the latter is pathological (see Sec. 2.5). 
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Scalar 


;f-ab 
'a 


q-)ab± 


^i(k)Ci 


-i9i(k)Ci7'^'' + ?9'i(k)kCin'^u^ + 2i?i(k)K«'' 


M^)r' 


^2(k)C2 


-i92(k)C27"'' + ??^(k)kC2u'^u^ 


??2(k)k''k* 



As one can see, there is a total contribution to the pressure and energy densities defined in (2.44). 
Summing up the two different contributions to the dipole moment, one concludes that for co-dimension- 
1 surfaces, it is necessary to add to the action composed of the terms (2.40) and (2.55) a 1st order 
contribution of the form 

^{i)[X^\= I V^(M^)r' + M^)^''^'')Kab= [ V^VfKab . (2.62) 

Here, T>°^ itself measures the most general response of a stationary fluid brane due to a 1st order bending. 
To second order there is also another set of contributions that can be added to the action which mix the 
fluid and elastic behavior^^. These are 



^3(k)k'^k^k^VaKb, , i94(k)k''VaK^ , T?5(k)k'^Vair , 



(2.63) 



aci,fe D± 



i?6(k)k'^i^^Vbk , ^7(k)7"'k''i?^,fe, 



However, the term 'i?6(k) can be related to the term 'i?3(k) by on-shell equivalences as in Eq. (2.57), while 
the last contribution '!?7(k) can be related to two other terms as we will explain below. It is then only 
sufficient to consider the first 3 contributions in (2.63). For this we need new machinery. For simplicity, 
we examine an action of the form: 



JWp+i 



(2.64) 



To study terms proportional to VaK^c it is necessary to define a new moment D^ of quadrupole nature: 



V 



,abc 



1 



6C 



(2.65) 



V-7 6VaKbc^ 

Using App. A for the correct variations we find the modified intrinsic and extrinsic dynamics of the 
fluid-elastic system: 

V.T'^^ - ZVaV^^^diKK - 2P'^'^^VdVaK^ = 2V''^'Wd\Kec + VJ)^''''R^cd , (2-66) 



T'^^'Kab + (3V,P°'"^ + 2V,P'"^^j KmK'', - 2VdV'''''KabK\ 



The extrinsic twist potential does not play a role in co-dimension- 1 surfaces since it vanishes. 



(2.67) 
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The boundary conditions are summarized in App. B. Inspecting Eqs. (2.66)-(2.67) it is easy realized that 
it does not fall into the class described by the equations of motion (2.34)-(2.35). Indeed, for this type of 
corrections to be properly accounted for, it would be necessary to understand quadrupole corrections to 
brane actions in the spirit of Sec. 3. Below, we summarize their contributions to the stress-energy tensor 
and corresponding quadrupole moment: 



Scalar 


'a 


'Tjabcl. 
a 


M^)C3 


-M^)C3r'' + ^3i^)^C3U''u' 


l?3(k)k'^k''k'^ 


^4(k)C4 


-M^)CaY^ + i9^(k)kC4n"n* + 2T?4(k)k^V('^K^)e 


??4(k)k(V^" 


^5(k)C5 


-^5(k)C57"' + ^'5(k)kC5n'^u'' + 2i?5(k)k^V,i^"^ 


^5(k)k'^7^^ 



The sum of the 3 contributions to the quadrupole moment can be written as 

Vf"" = t?3(k)k''k^k'= + ??4(k)k(V)" + T?5(k)k"7*'^ , 



(2.68) 



and encodes the most general response to corrections proportional to VaKf^c for co-dimension- 1 surfaces. 
Therefore, for co-dimension- 1 surfaces, one needs to add to the action (2.26), a 2nd order contribution 
of the form 



JWp+i 



(2.69) 



We end this section by noting that the terms 'i?4(k), ^^^(k) and i?7(k) are coupled to each other. To see 
this we contract (2.32) with 7°'^k in order to obtain: 



^^'k'R^.ba = k" (VbK - VaKb") 



(2.70) 



This is why we did not need to study the term 'i?7(k) separately. To summarize, the family of curved 
fluid membranes embedded in a curved background is parametrized by a total of 15 response coefficients, 
13 of which are independent. 

3 Multipole expansion as derivative corrections 

In this section we establish a relation between the action formalism of Sec. 2.2 and a multipole expansion 
of the stress-energy tensor for curved branes. This expansion, to pole-dipole order, is sufficient to capture 
all the corrections studied in Sees. 2.4-2.6 but not enough to capture the quadrupole corrections studied 
in Sec. 2.7. Here, we focus on the pole-dipole order and leave the extension to pole-quadrupole order 
for future work. The relation between this expansion and the corrections of the previous sections is 
particularly useful as it gives physical meaning to the corrections themselves and it allows to establish a 
precise link between gravity and the effective description of fluid branes. 
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3.1 Equations of motion 

The multipole expansion consists in expanding the stress-energy tensor of a curved brane in a Dirac-delta 
series, in the same spirit as in electromagnetism the electric current is expanded in order to describe finite 
thickness dielectric effects [43]: 

The stress-energy tensor (3.1) is symmetric, since it should be coupled to gravity, which implies B^^'^ = 
B^^^' as well as B'^'^'^ = B^^'^'i^ . The structure B^'^ is a monopole source of stress-energy while the struc- 
ture B^^'^P is a dipole source and encodes the finite thickness effects. Introducing higher order structures 
in the expansion (3.1) results in higher order moments such as the quadrupole moment. Further, to each 
structure in the expansion we associate an order parameter e such that B^^'^ = O (o) and B^^^ = O {e). 
More importantly, note that the expansion parameter e is not the same as the expansion parameter e 
introduced in Sec. 2 with the purpose of keeping track of the derivative expansion. For example, in the 
case of the hydrodynamic expansion for unbent black branes [12, 13, 14, 15], one can iteratively correct 
the metric to arbitrary order in e while still being at order O (o) in e. However, if the deformations 
of black branes involve strains along transverse directions to the worldvolume the corrections will be 
encoded in B^^^p [5, 16, 6], as well as in higher order structures, inducing a multipole expansion in e. In 
particular, corrections proportional to one copy of the extrinsic curvature are encoded in B^''^p . Note, 
however, that if B^^'^P contains a one-derivative terin, then the worldvolume effective theory that it gives 
rise to is of order O (e^) . This means, for example, that deformations of black branes can acquire dipole 
moments B^'^'' to order O (e) without affecting the effective theory to that same order. This is due to 
the fact that the covariant derivative in (3.1) is acting on B^^^p. This digression will become clearer as 
we progress in this section. 

The equations of motion for an object characterized by a stress-energy tensor of the form (3.1) can 
be obtained at the probe level, and in the absence of other external background fields, by solving the 
conservation equation 

V^f^"" = , (3.2) 

using the methods described in [43]. To express these equations in a nice form it is useful to decompose 
the structure B^'^p into tangential and orthogonal constituents to the worldvolume BP'^p = Iu^p^B'^^ + 
uP'au'^bB'^ with Bj^'^'"' = B^"" ' = 0^^ and to introduce the variables^^ 



The symbol _L here means that the space-time indices of the tensor are transverse to the worldvolume. For example, 
7\Bf *> = 0. 

^^Here we have gauged away the parallel components it''a_B'"'" since by the 'extra symmetry 1' [43] they can be gauged 
away everywhere on the worldvolume. On the boundary, these components may be non-zero but if we assume the absence of 
extra sources on the boundary they vanish. In any case, this extra structure can be trivially incorporated into our description 
by redefining B" — >■ _B° = _B° — u^ pU xWcB'' "^ and adding an extra term to the boundary conditions [43]. 
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In terms of these, the equations of motion can be written as [43] 

±^±%Va5^^" = , (3.4) 

Vb (m'^'ua^' - 2u\VaS^''"' + U^'cU%u\VaSP'"') = S^P'R^'.Xp ■ (3.5) 



Eq. (3.4) expresses the conservation of transverse angular momenta. In the point particle case, these 
equations reduce to the ones derived by Papapetrou for spinning point particles [63]. Eqs. (3.4)-(3.5) 
must be supplemented by a similar decomposition of B^" as B^'^ = B^ + 2M''^(,i?^ + u^ a^'^ bB""^ , where 
the first two sets of coefficients are related to the dipole structure via the relations [43] 

Sf = ±'^A^%VaiV^^'^ , B''^ = u''x^''p^b{s^P^ + N^P^) . (3.6) 



These relations imply that the only free coefficients in the theory are B""^ and B'^'^p. Eq. (3.5) can be 
rewritten in a simpler form. In order to do so, we decompose the term inside the parenthesis on the 
l.h.s. into tangential and orthogonal parts to the worldvolume and define a linear momentum tensor of 
the form^^ 

pup ^ / a6 _ u'^pU^VcSP^") U^aU\ + 1l\\-^ \V cS^^^ ■ (3.7) 



The linear momentum (3.7) is not necessarily symmetric. Further, it is tangential in its first index but 
not in its second, that is, l.PjyV'^^ = 0. Using this definition, the equation of motion (3.5) can be recast 
d la Carter^'^ [37]: 



7\va 



VxVP = S^'''R^cXp , (3.8) 



where the term on the r.h.s. can be seen as a force J-^ = S '^^R^cXp acting on the worldvolume due to the 
coupling to the background Riemann tensor. Moreover, note that when the Riemann curvature tensor 
vanishes, the linear momentum V^" is conserved along the surface. The set of Eqs. (3.8) can be split into 

h 

with n*^, leading to. 



two sets by projecting along the worldvolume directions with v!' n and orthogonally to the worldvolume 



VaV'''>-V'''K\i = S^P'R\xp , (3.9) 

V'^Kab' + yaV'+ujJjV''' =S^P''R\xp ■ (3.10) 

Here we have used the definitions V""^ = u°'yU^ ^jP^p and "P^* = u^yTi^ ^T"^^. These equations are written 
in the same fashion as in the work of Guven et al. [61], however, we have generalized it for pole-dipole 
branes and for curved backgrounds. Now note that in the case of absence of dipole effects S p^ = "P"* = , 
Eqs. (3.9)-(3.10) reduce to those obtained to 0th order in the expansion (2.7)-(2.8) upon the identification 

^^The reason for attributing the name of 'linear momentum' to (3.7) is explained in App. A. 
^^This sentence should be pronounced with French accent. 
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^ab _ rpab ^ rjij^g space-time stress-energy tensor that gives rise to the worldvolume theory (2.7)-(2.8) is 
the one given in (3.1) with B^'^P = 0. This is what is meant by the fluid being confined to an infinitely 
thin surface since the stress-energy tensor T'^'^ is localized there to order O (o). 

Before we understand the relation between the multipole expansion and the action variations of the 
previous section, it is important to provide the boundary conditions that arise from solving Eq. (3.2). In 
terms of the linear momentum (3.7) these can be written as 

(3.11) 

2Va (S'^'^'^r^avt) - VbV"'] lew.+i = , 



where we have used the definition of the boundary vectors introduced in Sec. 2.2. 

3.2 Relation to the action principle 

To connect the results of the previous section with those of Sec. 2 it is necessary to clarify the physical 
meaning of the different components of the structure B^^^p. This has been done in [5] and we review it 
here. If we focus on flat space-time and on uniform p-branes extended along all x", ..., x^ directions we 
can evaluate the total angular momentum on the transverse plane labelled by the indices //, v as: 

where S is a spatial slice of the background space-time. This leads to the introduction of a spin current 
ja^u g^(.j^ that j"^*^ = 2B^°' . Further, we can evaluate the worldvolume dipole moment of the brane as 

j^abp ^ I d^~^xfp''u^''uJ'xP = [ dV^p^B'^'' , (3.13) 

leading to the introduction of a worldvolume dipole density d°'^P = B'^ . This last one, as we will see 
below, can be interpreted as the bending moment of the brane. Given these definitions we can rewrite 
the tensors 5^'''' and NP""" introduced in (3.3) as 

gl.ua ^ }_■a^,u _ ^abi^^^u]^ ^ ^(.ua ^ ^ab{^,^u)^ _ (3^^4^ 

We will now focus on different aspects that relate the multipole expansion (3.1) with the action (2.26). 

The dipole moment 

Here we consider the case where the transverse angular momenta (3.12) vanishes, that is, j"'P'^ = 0. In 
this case, the linear momentum (3.7) can be written as 

p./. = (f°^ - d'^'^K^^ uPau\ + uWPxVad''''^ , (3.15) 
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where we have defined the symmetric tensor T""^ = 3°"^ + 2S-'^'^^ K^' d. With this, the equations of motion 
(3.9)-(3.10) can be written as^*^ 

Vaf^^ = -u\Vaycd'"'^ + d'''''R\ci , (3.16) 

f''''Kj = -n\VaVbd^'"'-d'''^R\jb ■ (3.17) 

Comparing this with Eqs. (2. 34), (2. 35) when the spin current vanishes, that is S^ij = 0, we find that 
both are equivalent upon identifying 

which in turn imphes B""^ = T""^ — 2d^°''^^K^'ci- The boundary conditions (3.11) and (2.31) can also be 
seen to be equivalent upon the same identification. With this we have shown that all corrections to brane 
actions quadratic in the extrinsic curvature can be accounted for by a multipole expansion to pole-dipole 
order of the stress-energy tensor. Moreover, the conservation equation (3.4) reduces to the integrability 
condition 

(f^^'^Kab^^ = . (3.19) 

This condition is automatically satisfied for the actions we consider due to the form of T)""^^ given in 
(2.45). Indeed, the form of (2.45) is what is expected for the bending moment of thin membranes 
[57, 5], an interpretation which is now justified due to Eq. (3.13). In fact, when a rod is bent, a varying 
concentration of matter across the transverse directions induces a bending moment which is proportional 
to the Lagrangian strain (2.14) [57, 5]. Further, note that this analogy with classical elasticity is direct 
when one deals with co-dimension- 1 surfaces for which the transverse index in the extrinsic curvature 
can be omitted [57]. 

The spin current 

Now we consider the case for which the worldvolume dipole moment vanishes d = but the transverse 
momenta is non-zero. In this case, the linear momentum (3.7) takes the form 

V''^' = B''\^'au''b + u\l.%,fP'K\, . (3.20) 



In turn, the intrinsic equation of motion (3.9) reads 

1 
2 

while the extrinsic equation (3.10) takes the form [5] 



Vai?"' = --J^O,"*^- , (3.21) 



B-'Kab' = -n\Vb (ja'^^K'^'j) + lr^'R\kj . (3.22) 



These equations have been previously obtained in [5] but in a different form. 
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Comparison of this last equation with Eq. (2.35) when dipole effects are absent P"** = leads us to 
identify 

^ab ^ j,ab ^ -aij ^ _25«U . (3.23) 

This identification is sufficient for the intrinsic dynamics (2.34) and boundary conditions (2.31) to match 
those given by Eq. (3.21) and Eq. (3.11) when Eq. (3.4) is imposed. The conservation equation (3.4) 
reduces instead to the conservation equation of the spin current 

i±^±%V«r^^ = . (3.24) 

This is the reason why we have discarded the term tui(k) in Sec. 2.5 as it does not satisfy this equation. 
Given this identification we have shown that corrections quadratic in the extrinsic twist potential^^ can 
be accounted for by a multipole expansion of the stress-energy tensor and that these are interpreted 
as the fluid-elastic system acquiring motion in transverse directions to the worldvolume by means of 
Eq. (3.12). Both the integrability condition (3.19) and the spin conservation equation (3.24) may be 
derived as constraint equations from the action (2.26) in the spirit of [61]. We leave this exercise for 
when we analyze in detail spinning corrections to black holes [60]. 

The 'extra symmetry 2' 

The stress-energy tensor (3.1), truncated to order e, enjoys of a perturbative symmetry coined the 'extra 
symmetry 2' by the authors of [43] which results from the freedom of shifting the worldvolume surface 
by a small amount, that is, 

Xf^ia") -^ Xf^ia") = Xf^ia") + e^(a'^) , (3.25) 

where e'^ is an infinitesimal shift vector of order i. The interpretation of this symmetry can be understood 
if one remembers that the expansion (3.1) delocalizes the stress-energy tensor by giving the surface a 
finite thickness. Since the thickness is finite, the worldvolume can be placed anywhere inside the surface 
[43] . This symmetry acts on the structures B^'^ and B'^'^p as^^ 

52B^"' = -B^"'u%Vaef - 2B^^^'T''\peP , 

(3.26) 
d2B''^P = -B'"'eP , 

and leaves the stress-energy tensor (3.1) invariant to order i. Along worldvolume directions, this transfor- 
mation coincides with worlvolume reparametrizations^^ but it is non-trivial along transverse directions. 



^^This also holds when the corrections are proportional to only one copy of the extrinsic twist potential. 

^"in order to obtain this transformation rule from (3.1) one should use (2.3) and the fact that dpB'^"''' — since B^"'' is 
a function of cr" and not of the space-time coordinates. Further, note that any scalar or tensor which is a function of the 
space-time coordinates does not transform under a shift of the worldvolume surface. See Ref. [43] . 

^^This is true except at the boundary when extra boundary sources are present [43]. 



25 



To see this directly in the equations of motion we separate the transformation into parallel and orthogonal 
parts such that e^ = e^n^i + u^a^""- Then, according to (3.26) we have that along orthogonal directions 



-'2 



m 



«6 = - (u%m''^ + n(%m'')^) V^ (e^n^^i) , d^cP^^ = -m"¥ , d^f'^ = O [e^) . (3.27) 



Further note that due to the identifications (3.18) and (3.23) the stress-energy tensor T°-^ transforms as 
(52 T"^ _ g^YYiO-b _ rn'^y°- K^> cie'^ . Since the transformation oi j°-'^^ is of higher order it suffices to look at the 
equation of motion (3.17). Using the transformation rule (3.27) together with (2.18)^^ it is straightforward 
to check that the equation of motion (3.17) is invariant under the 'extra symmetry 2'. This can also 
be seen at the level of the action (2.47). In this case, note that m°'^ can be replaced by T?^s^ due to 
the identification (3.18) and (2.41). Therefore the variation of the dipole moment can be written as 
^J_ja6i _ _ j'afegi^ 'pj^jg nieans that, by virtue of the definition of the dipole moment (2.27), a term of the 
form —T^^.eiKab must be added to the action. Now, from the form of the action (2.47) it is necessary 
to know how — -^— 7'P(k) transforms. Using Eq. (2.3) this is simply 

5i{-V^P{\^))=T^^f^Ka,' . (3.28) 

Clearly, due to the transformation rule (3.28) together with the transformation of 5^d"** the action (2.26) 
is invariant under the 'extra symmetry 2'. Now consider the case in which the expansion parameter e 
is controlled by the radius of curvature of the worldvolume. In this case we can write e*((j'*) = k(a"')K^ 
for some arbitrary function k{a°'), leading to (5^(i"''* = —T?'^Jk{a°')^'^'^Kcd. This means that the Young 
modulus (2.46) picks up a gauge dependent part of the form 

abed u( ^CL\ I rpab ^fid , r-ncd ^,ab 



ya,ca ^ _^(^a) iT?'^''' + T^'^'^a" 

) ' , (3.29) 

= -fe(a") ( -2Ao(k)7'^S"^ + A'o(k)k(n"nV^ + n^A"^) 



where we have used the explicit form of T?z. given in Eq. (2.10). This indicates that the elastic con- 
tributions Ai(k) and A4(k) given in (2.40) pick up gauge dependent terms. Similarly, for surfaces of 
non-vanishing mean extrinsic curvature, this can be understood as the pressure i-*(k) picking up a de- 
pendence quadratic in the extrinsic curvature such that, 

P(k) ^ P(k) = P(k) - T^^^-f'''k{anKab'K,d^ . (3.30) 

Thus, a natural choice of gauge is one for which the pressure P only depends on k^^. It is worth noting 
that since the variation of j"*-' is of higher order, the equation of motion (3.22) is not invariant under 
the 'extra symmetry 2'. This in fact means that the truncation d°^* = of the pole-dipole equations 
(3.4)-(3.5) is not gauge invariant. To make it gauge invariant one should add an extra d term of the 



^^When using Eq. (2.18), $' should be replaced by e\ 

■^^Note that the inclusion of the hydrodynamic contributions (2.55) does not spoil the invariance of the action (2.26) under 
the 'extra symmetry 2'. This is because since these scalars only contribute to the monopole stress-energy tensor, their 
variation is of order O (e^)- 
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form d"** = — ?7T,"*e*. Relating this with the action (2.53) which accounts for spin corrections, it means 
that one should in fact add a term of the form {l/2)y°'^'^'^Kab^Kcdi in case of non-vanishing mean extrinsic 
curvature or a term d Kabi generically. However, we can always choose a gauge for which P is only 
dependent on k, which is the gauge choice leading to (3.22), or equivalently, choose a gauge for which 
d"^* = 0. 

3.3 Construction of conserved charges 

In this section we show how to construct conserved currents and charges for systems obeying the pole- 
dipole equations of motion (3.4)-(3.5). We begin by using the technique developed by Carter for 
geodynamic-type branes [36]. This consists in finding a conserved surface current V^ , by definition 
purely tangential J-^i^V^ = 0, such that 

7\VaP£ = . (3.31) 

The surface charges constructed from such current would then be conserved charges of the system. We 
will now show that the same type of ansatz as the one used by Carter can be extended to a large class 
of pole-dipole branes. We take the surface current to be of the form 

Va = Vi^k^ + J^^i^PV f,kp , (3.32) 

for an arbitrary space-time Killing vector field k^ and with J-^^Ti'^'^P = 0. Here we have used the definition 
of the linear momentum (3.7). We now introduce the ansatz (3.32) into Eq. (3.31) and use Eq. (3.8), 
leading to 

k^ (^S^P%\ + S'^^^) i?^,Ap + [V^"''^ + VeS^l'^^]) V.k^ = , (3.33) 

where we have used the fact that for any Killing vector field Vj^V^kp = Rp^,y\h . Looking at the term 
proportional to the Riemann tensor, it seems that the correct choice of Ti'^'^P is 

E'^^'P = -S^'f'%\ , (3.34) 

as long as 'pl'^w = —S/c'^'^^'^'^' is satisfied. We now focus on the class of pole-dipole branes for which the 
integrability condition (3.19) is satisfied. In this case, a simple exercise using (3.7) allows one to show 
that 

pM = Vc^t^'^l^ - ^^A^'^pVa^^^" = Vc^I'^^l^ , (3.35) 

where we have used Eq. (3.4) to eliminate the second term on the r.h.s. above^''. Hence, the choice 
(3.34) leads to a conserved surface current of the form (3.32). For the case derived by Carter in [36], 



■^^It may be useful to write the l.h.s. of Eq. (3.4) explicitly in terms of the variables introduced in (3.14). This is simply 
the sum of Eq. (3.24) with Eq. (3.19). 
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where j^^'^ = 0, we have that S'^'^^ = —(F^p. Furthermore, note that the integrabihty condition (3.19) is 
satisfied for all the actions we considered in Sec. 2. 

Carter suggested in [36] that charges obtained from the relevant surface integral of the current (3.32) 
are the conserved charges of the pole-dipole brane. The naive way to implement this would be similar to 
the 0th order case (2.23), that is, a generic charge Qk associated with a background Killing vector field 
is simply given by the integral of the surface current (3.32) over spatial slices of the worldvolume 

IQkl = / dV(^)V!^n, . (3.36) 

JBp 

However, we can test if this is indeed the case, since we have the full space-time stress-energy tensor 
(3.1). Since (3.1) is symmetric, any current of the form P^ = T^^\<i^{x'^) will be conserved in the full 
space-time, that is, 

V^P^ = . (3.37) 

Under the assumptions that the timelike Killing vector field ^^ is hypersurface orthogonal with respect 
to the space-time metric g^y and that ^^ is parallel to the worldvolume timelike Killing vector field E,"", 
i.e. ^^^ = u^ai"" 1 which is also assumed to be hypersurface orthogonal with respect to the worldvolume 
metric, we can write the total conserved charge as^^ 



IQkl 



dV^^)B^''n^K+ i dV^p)RoB^^''f>V,[^ 



(3.38) 



where all quantities should be evaluated at the worldvolume surface x° = X°(o""). Under the assumptions 
mentioned above we can write the normal vector Up to a spatial slice as rip = —Rq6p^, where t denotes 
the time coordinate of the background space-time. Therefore, the difference between the charges (3.36) 
and (3.38) is given by 

IQkl - IQkl = 2 / dV(^p)Bf'''^'npKdplnRo . (3.39) 

JBp 

A few comments are now in place. In general, the charges computed from conserved surfaces currents 
differ from the charges computed from conserved space-time currents. In flat space we have Rq = 1 and 
hence the charges agree. Clearly, from the form of (3.39) the surface charge computed from (3.36) is not 
taking into account the variation of the redshift factor Rq along orthogonal directions to the worldvolume. 
In the 0th order case, for which the brane is infinitely thin, this is not a concern but when the brane 
acquires a finite thickness it is natural to take this effect into account. Given (3.39) one should either 
modify the definition of the current (3.32) such that 

^k^^k=^k-2 5''^%9plni2o , (3.40) 



■^^Note that the expression for the conserved charges Eq. (3.38) differs from the one obtained in [5]. The expression 
obtained in [5] is only correct in flat space as it does exhibit the full invariance under the 'extra symmetry 2'. This issue is 
clarified in App. A. 
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at the expense that V^ no longer satisfies the surface current conservation equation (3.31), or simply 
subtract the additional term in (3.39) to the definition of the charge (3.36) itself. Note that the definition 
(3.40) implies that the current "P^ does not satisfy in general the tangentiality condition, i.e., J-^^V^ / , 
suggesting that this is a finite thickness effect. It is worth noting that since (3.1) is invariant under the 
transformation rule (3.26), so are the charges (3.38) computed from it for any choice of k. In order to 
see this explicitly one should use (3.26) in (3.38) together with (2.3)^^. On the other hand, the charges 
computed from the naive method (3.36) are not invariant under the 'extra symmetry 2' in general, only 
in the case of flat space. It is worth noting that as in (2.23), the choice of Killing vector field k in 
Eq. (3.38) results in either the total mass, angular momentum along worldvolume directions or angular 
momentum along transverse directions to the worldvolume of the fiuid-elastic system. Moreover, in the 
limit where both dipole and spin effects are turned off, the charges computed from Eq. (3.38) agree with 
those computed from (2.23). As a final comment, we note that the conserved surface current (3.32) can 
be obtained directly from (2.26) by requiring the action to be invariant under space-time translations 
along Killing directions. We show how this is done for the case j"'^'^ = in App. A. 

4 Matching with gravity 

In this section we apply some of the results of the previous sections to the bending of neutral black 
branes in pure Einstein gravity. It is well known that the stress-energy tensor of black p-branes takes 
the form of a perfect fiuid [19] and if the brane is subject to a long-wavelength perturbation then the 
stress-energy tensor (2.10) receives corrections at each order in the derivative expansion. If the type 
of perturbation is along the worldvolume directions then the corrections are generically of dissipative 
nature and do not fall into the class studied in Sec. 2. However, if the perturbations are along transverse 
directions and lead to stationary configurations then the system (2.26) could potentially describe the 
effective behavior of the perturbed brane. Indeed, this is the case to 0th order (2.2) as demonstrated in 
[17, 16]. The interest in this type of stationary perturbations comes from the several applications of the 
blackfold method [17, 18, 19] to the construction of higher-dimensional black hole solutions. The method 
consists in taking the metric of a boosted black p-brane and wrapping it over a submanifold Wp+i of 
characteristic curvature R. The end product of such methodology is the perturbative construction of 
higher-dimensional black holes for which their near horizon geometry is that of a bent boosted black 
p-brane to a certain order e = ro/R in the perturbative expansion. Here tq is the thickness of the brane, 
which for Schwarzschild branes coincides with the horizon radius. To be more precise, it is instructive to 
write the metric of a boosted Schwarzschild black p-brane mD = n + p + 3 space-time dimensions 

d^o) = ( 7a5(X^(a")) + J>A-i^,(a")n,(a") ) da^da" + ^.^^^^^ + r'cinl^,^ + (4.1) 

^ 'J- ~.r? 



■^^In fact, any scalar functional of the form T[f] = J d'^'a;-^— gT'"^/,i^(a-°), where f;_w{x°') is an arbitrary tensor field of 
compact support, is invariant under the transformation (3.26) [43]. See App. A to understand how (3.38) is derived. 



29 



Here we have promoted the horizon radius tq as well as the boost velocities u" and brane worldvolume 
metric 'jab to slowly varying functions of a"" over the submanifold Wp+i. When promoting the various 
fields to functions of the worldvolume coordinates, the metric (4.1) is in general no longer a solution of 
Einstein equations and should be corrected by including terms proportional to the derivatives of ro, u°' 
and 7af,- To 0th order in e the metric of the black hole solution constructed from wrapping Scharzschild 
p-branes is given by (4.1) and since, as it will be explained below, its stress-energy tensor is of the perfect 
fluid form, the effective dynamics are described by a system of the form (2.2) with equations of motion 
(2.7)-(2.8). In fact, this has been shown to be the case directly from Einstein equations [17, 16]. Here 
we are interested in perturbations which are first order in derivatives of the metric "fab along transverse 
directions. In this case, to first order in e, a small perturbation h^^j should be added to the metric (4.1). 
Generically, it can be put into the following form [16] 

„2 _/„ ot- i „a , 'O.. „. \ ^^aib , "' , „2,/)2 , „2 „■ 2 /j ,r.2 



dsl^) = [riab- "^KabT COS e + -^UaUb da" da" + -^ + r^dO^ + r' sin^ Qd^^. 

^' \ r^ J 1 _ ^ '^ > (4.2) 

+ V(r, e)dx>'dx'' + O {r^/R^) . 

Here we have used the label i to indicate that the perturbation is taken along a single transverse direction 
to the worldvolume and also introduced rjab, which is the flat metric on the worldvolume. As explained 
in [64, 16], perturbations in each individual transverse direction decouple from each other at this order. 
The important piece of knowledge about the perturbation hfj^u is that it is a dipole type perturbation 
hfj,u {r,0) = cos Oh^u (r ) and that all of its components are proportional to the extrinsic curvature Kab [16] . 
For the metric (4.2) to be regular it must satisfy the equations of motion (2.7)-(2.8), which can be shown 
through the usage of the method of matched asymptotic expansion [17, 16]. The metric (4.2) is valid 
to order 0(e) and promoting all fields to functions of a"" as in (4.1) allows for iteratively correcting the 
metric of black holes constructed from wrapping Schwarzschild p-branes. However, as argued in Sec. 2.3 
and shown in [17, 16], there are no corrections to the asymptotic charges, angular velocities, entropy and 
temperature to order 0{£) from (4.2). However, as explained in the beginning of Sec. 3.1, introducing 
corrections proportional to the extrinsic curvature induces a multipole expansion of the stress-energy 
tensor (3.1). Therefore, even though there are no corrections to the charges, the metric acquires a 
bending moment [5, 16, 6] to this order, which we will analyze below. 

There is an exception to the absence of corrections to the charges to this order, which is the case n = 1, 
related to the presence of backreaction effects [17, 16]. As argued in [5], there are two types of corrections 
that (4.1) can be subject to: backreaction corrections and curvature corrections. Newtonian estimates 
[5] indicate that curvature corrections become more important when n > 2. Curvature corrections (or 
elastic) are included in the formalism of Sec. 2 since the resulting equations of motion satisfy stress-energy 
conservation (3.2), while backreaction effects are not. In the cases where backreaction is subleading, the 
perturbation (4.2) is said to be a pure bending. Therefore, it is only for the cases n > 2 that one should 
expect an effective description of the form (2.26). The procedure of iteratively correcting the metric by 
introducing strains along transverse directions to the worldvolume has only been completed to order O (e) 
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leading to the metric (4.2) and, due to the current state of affairs, there is no data available to O (e^) in 
this long-wavelength perturbation. Nevertheless, it is possible to obtain information about (2.26) from 
the metrics (4.1) and (4.2) as we will describe below. 

4.1 0th order metric: thermodynamic fluid variables 

To 0th order the metric describing the fluid-elastic system is the one presented in (4.1). The stress- 
energy tensor obtained from the ADM formalism or the Brown- York prescription is of the perfect fluid 
form (2.10) where the pressure and energy density take the following form [18] 

The local fluid temperature T and local entropy density s can be obtained from (4.1) by reading off the 
surface gravity and the horizon area respectively, 



n 



n 



("+!) ^n+1 



The Gibbs-Duhem relation e + P = Ts and the first law of thermodynamics de = Tds are obeyed 
due to the form of (4.3)-(4.4). Due to the relation between global temperature and local temperature 
T = k?" mentioned above Eq. (2.12), one can find the relation between the thickness tq and k such that 
(47rr)ro = nk and hence the dependence of the pressure P given in (4.3) on k , 

-«-^(i^)"^"^ (-) 

The dependence of P on k is in fact enough to predict the energy density and entropy through Eq. (2.10) 
and Eq. (2.12) . Using the identification (2.11) one has that Ao(k) = — P(k), and hence that AQ(k) = 
— nP(k)/k. Therefore, using (2.10) we find 

T^^^ = P(k)7'^^ - nP{\L)u''u^ , (4.6) 

which is the stress-energy tensor for the fluid (4.3)-(4.4) with equation of state e = —{n + 1)P. To 0th 
order, we have full predictability of black hole masses, angular momenta and entropy using (2.23) and 
(2.24) since we know all the local microscopic properties of the fluid (4.3)-(4.4) . This has been used to 
show the existence of several new black hole solutions [65, 66, 54, 55]. 

4.2 1st order metric: Young modulus 

As mentioned above, there are no corrections to the asymptotic charges and local thermodynamic poten- 
tials to order O {e) of the metric (4.2) when n > 1, however, the bending moment D"''* given in (2.45) is a 
0(e) correction since it is only proportional to one copy of Kab ■ The bending moment can be measured 
from the metric (4.2) using the methods of [5] and briefly, consists in finding the stress-energy tensor of 
the form (3.1) that sources the metric (4.2). This was first done in [5] for neutral black strings bent into 
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a circle and later generalized for Schwarzschild black p-branes bent into an arbitrary shape in [16]. The 
bending moment found in [16] is of the form (2.45) as expected from classical elasticity theory with a 
Young modulus given by 



(4.7) 



yabcd =p(k)r2(k)e(n) ('^^7«(^7'^)^ + 2n('^7'')(%'^) + ^^^^-^u^u^^uA 
- k P(k) e(n) (27"^^"^ - ri iu^'u^-i''^ + n'^u V^) ) , 
where A; is a constant and the function ^{n) takes the form 

_ ntanWn)r(H±i)^ 

' r(H±2)= 

Note that the function ^{n) evaluates to zero when n = 1 and diverges when n = 2, in agreement with 
the expectation that elastic corrections are subleading when compared to bakreaction corrections for the 
cases n = 1,2. One of the key results in this work is the prediction of the structure of (4.7) from (2.46). 
Indeed comparison of (4.7) with (2.46) using the results of Sec. 3.2 leads to the identification 

A.(k)^-.P(k),(n) , A.(k)^ WIM , A3(k)^MM^ 



(4.9) 



fcnP(k)g(n) 3n + 4 P(k) r2(k) ^(n) 

MW = r-2 ' -^sikj — 



2(n + 2) 



A few remarks are now in place. Note that the elastic contributions Ai(k) and A4(k) are gauge dependent. 
Inspecting the second line of (4.7), we see that it has the same structure as the gauge dependent part of 
the Young modulus defined in (3.29). This comparison allows us to identify k{a°') = k^{n). Furthermore, 
it is not surprising that the scalars associated with Ai(k) and A4(k) are only gauge dependent. This is 
because since the metric (4.2) is valid only to order O (e), it cannot probe the hydrodynamic contributions 
(2.55) since these are second order corrections to T'^^ . Therefore, since TZabcd can be neglected to this 
order and since the measurement of (4.7) was done in flat space, according to Eq. (2.58) we have that, 
£1 = £2 and £3 = £4. It is worth noting that measuring (3.1) from the metric (4.2) does not uncover 
any of the corrections to t or 11 given in Sec. 2.4 and Sec. 2.6 since these are of order O (£^). In 
fact, besides the measurement of P'^**, from (4.2) one also obtains B""^ as defined in (3.3) which takes 
the same form as in (4.6). In other words, to order 0(e), the perturbation (4.2) does not affect the 
monopole contribution to the stress-energy tensor (3.1). This means, for example, that one cannot see 
the full invariance of (3.1) under the 'extra symmetry 2' to order O (e). 

4.3 Elastic corrections to black rings 

As mentioned towards the end of Sec. 4, the procedure of iteratively correcting the metric (4.1) in a 
derivative expansion has not been completed to order O (e^). However, the measurement of P"''* is all 

32 



that is necessary to predict the structure of r"^ using the table given in Sec. 2.4. One may assume that if 
the deformation to the metric (4.1) is a pure bending then the only excited modes are those presented in 
(2.40) and the hydrodynamic contributions ^"2(1^) and U3(k)^^. However, we can turn the hydrodynamic 
contributions off by considering the simple embedding of a string bent into a ring, for which TZabcd = 0, 
placed in a flat background. This means that the elastic contributions Ai(k) and A4(k) continue to only 
be gauge dependent. Further, due to the homogeneity and isotropy of the ring solution we take k to be 
constant along the worldvolume, which means that all the contributions from the hydrodynamic term 
t'i(k) will vanish (see table in Sec. 2.6). Given these facts, what we would like to test is if 1st order 
data is enough to predict black hole charges to 2nd order in the derivative expansion, at least for specific 
cases. For this purpose one would require the knowledge of the thermodynamic quantities (4.3) and, in 
particular, of the dependence of the pressure on k to order O (e^). From here on we will assume that 
this functional dependence to order O (e^) is the same as that given in (4.5). The effective action for this 
configuration would then just be that given by (2.47) with the identification (4.9). Using the definition 
of the modified pressure (2.44), this is simply 

F[X^] = - f V^V{k,Kj) . (4.10) 



n+l 



We parametrize the flat background as 

ds'^ = -dt^ + dr"^ + r^dc/.^ + ^ dx^ , (4.11) 

and place the ring at r = i? such that a^ = t = t , ,a^=(p and x* = 0. The induced metric is 
just 'jabdcr'^da = —dr"^ + R^d(p^. The only non- vanishing component of the extrinsic curvature is given 
by K^^^' = —R. This implies, due to the form of the linear momentum (3.15) that, "P"* = and the 
equations of motion (3.10) reduce to 

V'^Kab' = . (4.12) 

Since that there is only one non- vanishing component of the extrinsic curvature, (4.12) is actually just 
'P'P't' = 0. Instead of computing all the components of the stress-energy tensor involved in V'^'^ using the 
table given in Sec. 2.4, it is more practical to simply compute the scalar V (\i.,Kab^) and vary (4.10). In 
order to do so, we choose a gauge for which the leading order pressure P does not depend on the extrinsic 
curvature, which in turn implies that we choose /c = in (4.7) (see Sec. 3.2). Using (2.46) together with 
(4.9) we find 

V (k, AW) = (P(k) - ^ - A3(k)02 - X^{k)n'RA . (4.13) 

Varying now the action (4.10) leads to the equation of motion (4.12). This can be solved perturbatively 
in the manner Q = 17(0) + r2(2)e^ and T = T^g) + ^(2)^^! where e^ = r'^/R'^. Here 0(o) ^^^ ^(o) ^or 



■^^There is no Gauss-Codazzi-type equation tliat relates any of the elastic contributions (2.40) to the hydrodynamic 
contribution i;i(k) so if the deformation is pure bending this term is not, in principle, excited. 

33 



this configuration can be obtained from (2.2) together with (4.4) and are given by the relations ^(o)-R = 
l/\/n + 1' and 47r-v/rr+T'T(o) = n2 [65]. Introducing this decomposition in the equations of motion (4.12) 
leads to a solution to 0,i2) of the form 

Note that the end result (4.14) does not depend on the correction Tm). Given the corrected rotation 
velocity fi we can proceed and compute the total mass and angular momentum using (3.36). For this we 
need the conserved surface current (3.32), which for this particular case takes the form^^ 

Pl^ = p'^f^k^ - d^Mv ^kp . (4.15) 

Thus, the mass of the fluid-elastic system is given by 

M= f dV^p) V^^ = 2TrR T^^ = 27r i? P(o) ( (n + 2) - ^^"^ ^ ^^^ C{n) eA , (4.16) 

Where P(o) is the modulus of P(k) given in (4.5), when Q = il.(Q\. The angular momentum along the 
direction (p reads 

^ = -/^^^(p)(^"^xi^^ + Id^'-^adrxT) =2-^^P(0) V^^TT (l + ''^^,l\\''h in)e') . (4.17) 

As mentioned in Sec. 2.2, an entropy current formalism has not been developed for the system (2.26) so 
we do not have a first principle method to compute the total entropy of the fluid-elastic system. However, 
this is not necessary in order to check if 1st order data can account for black hole charges to 2nd order 
neither if the action (2.26) with the identification (4.9) is the correct effective description of higher- 
dimensional black rings in the blackfold regime. Since the action (4.10) must be interpreted as the free 
energy (see Sec. 2.1) we can obtain the product (TS)^ using (TS)^ = M — VlJ —F. Furthermore, since the 
charges M and J are supposed to be the ones associated with a stationary black hole in asymptotically 
flat space time, they must satisfy the Smarr relation 

(T5) =i!l±llM-f7J . (4.18) 

^ '^ {n + 2) ^ ' 

A simple exercise tells us that {TS)^ = {TS)2 + O (e^) for any value of n. Note that this result is again 
independent of the correction Tm) . Therefore, we conclude that 1st order data in the case of black rings is 
enough to have full predictability of black hole charges to 2nd order. We will discuss in the next section 
other possible cases for which 1st order data is enough to predict charges to 2nd order. 



■^^Note that in flat space-time, charges computed from the conserved surface current (3.32) are equivalent to those computed 
from the space-time current using Eq. (3.38) . 
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5 Discussion 

In this section we summarize the main results found in this work and discuss various open problems. We 
begin by stating the key results. In Sec. 2 we have found the most general action quadratic in the extrinsic 
curvature as well as in the extrinsic twist potential and in second order worldvolume derivatives. For 
co-dimension- 1 surfaces it was required, to the same order, to consider terms proportional to worldvolume 
derivatives of the extrinsic curvature. It was shown that the equations of motion obtained from this type of 
actions provide a relativistic generalization of classical elasticity theory of thin membranes when bending 
effects as well as spin effects are taken into account (2.36). Since the well studied case of fluid membranes 
is described by an action of the type (2.26), our work ended up generalizing the Helfrich-Canham bending 
energy (1.1) to the case in which the fluid living on the membrane is stationary and for non-trivial response 
coefficients. In such cases, for co-dimension- 1 surfaces, there exists an extra contribution to 1st order 
in the derivative expansion (2.61) and 6 extra contributions to second order (2. 40), (2. 63) than those 
considered previously in the literature. These response coefficients could potentially be measured in a 
physical experiment involving fluids moving on cellular membranes. 

In general, the results of Sec. 2 indicate that the study of hydrodynamics of fluids living on surfaces of 
arbitrary co-dimension is of increased complexity when compared to the hydrodynamics of space-filling 
fluids. In fact, for neutral stationary fluids it was found a set of 3 response coefficients [46] while for 
non-dissipative fluids a total of 5 transport coefficients [50]. Here instead, for surfaces of co-dimension 
greater than one and according to certain assumptions^^, we have found a total of 9 independent response 
coefficients, while for co-dimension- 1 surfaces we found a total of 13 independent response coefficients. 
In this counting we have ignored the spin contribution roi(k) as it does not fit into the pole-dipole 
approximation of Sec. 3, which suggests that it violates spin conservation. We have also shown in 
Sec. 2.6 how the elastic and hydrodynamic modes couple to each other using the Gauss-Codazzi equation 
and in which conditions we can regard each of the different contributions as independent. 

The work presented in Sec. 2 also indicates that all the techniques and methodologies used to study 
hydrodynamics of space-filling fluids can be applied to the case of fluid branes when the elastic modes 
(2.40) are taken into account. The requirement of stationarity can be relaxed and the method used here 
can also be applied to non-dissipative fluids as in [50]. In this case, the material space introduced in [50], 
where the fluid variables are defined, must be formulated with respect to the p spatial directions of the 
worldvolume Wp+i. 

We have found in Sec. 3 that all the corrections studied in Sec. 2 can be accounted for by the formalism 
of Vasilic-Vojinovic [43] where a multipole expansion of the stress-energy tensor is carried out to pole- 
dipole order, except for the terms (2.63) which require an extension of these ideas to pole-quadrupole 
order. Indeed, the formalism constructed by these authors provides the most general equations of motion 
that take into account finite thickness effects of curved branes, regardless of the existence of any underlying 
effective action. Having established this connection, we note that most of the finite thickness corrections 



We remind the reader that in Sec. 1 we stated that terms constructed with the Levi-Civita symbol would be ignored. 
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to brane effective actions considered in the literature [33, 34, 45, 36, 37, 38, 61] fit into the formalism 
of Sec. 3. In particular, we established a precise map between extrinsic curvature corrections and the 
bending moment of the brane as well as extrinsic twist corrections and the spin in transverse directions 
to the brane. It would be interesting to understand if this map can be useful for the effective description 
of long strings [67] . 

In this work we have generalized extrinsic curvature corrections to brane effective actions to the 
case of non-extremal branes and connected the formalism of Vasilic-Vojinovic [43] with the formalism 
of Carter [37] and the formalism of Capovilla-Guven [38]. The establishment of this connection is of 
particular usefulness as it allows us to connect effective theories of fluid branes with gravity. In Sec. 4 
we have shown that the Young modulus measured from bending neutral black branes falls into the class 
predicted by the effective action analysis of Sec. 2 and used this fact to predict the corrected horizon 
angular velocity for thin black rings. This gives further motivation for studying the elastic expansion 
of higher-dimensional black holes via the blackfold approach since now, uncovering the various response 
coefficients of stationary black branes simultaneously implies the uncovering of the possible response co- 
efficients and the structure of the free energy of real fluid membranes. This fact can be put into a broader 
context: this is another instance where the study of gravitational physics can shed light into the physics 
of an apparent unrelated system. Therefore, since thin black branes seem to behave like the membranes 
of living cells, gravity can be used as a laboratory for uncovering properties of fluid membranes. 

We now turn to several of the open problems encountered in this work: 

The spin/elastic contributions ■we missed: As mentioned in Sec. 1, we have ignored terms con- 
structed from the Levi-Civita symbol, in particular, we have ignored dimension-specific corrections of the 
form 

k^e^LOa^ , e''\jnab'' , e^\jRab'' , e^'i.KajKb'^ . (5.1) 

The second term above has been studied in the context of cosmic strings [37]. Further, the last 3 
contributions are coupled to each other due to the Ricci integrability condition (2.33). It would be 
interesting to understand if this last set of terms can be accounted for by the formalism of Sec. 3 to pole- 
dipole order. The first of these terms is of particular interest as it gives rise to the spin current measured 
from Myers- Perry branes [5]. Here, we have introduced the Levi-Civita symbol on a transverse two-plane 
labelled by the indices i,j. To make this statement precise, note that for this case 5"*-^ = k^e*-'. If one 
compares this with the result for Myers- Perry branes found in [5], this is exactly the type of correction 
needed for describing black holes spinning in transverse directions to the worldvolume. For the case of 
Myers-Perry branes we have access to the local thermodynamic quantities and hence in principle have 
full predictability of black hole charges when backreaction effects can be neglected. This problem will be 
addressed in a future publication [60]. 
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The entropy current and charged fluid branes: In this work we have lacked a first principle 
computation of the entropy for the fluid-elastic system (2.26). However, using the methods of [68, 46] 
together with some inspiration from the formulation of viscoelasticity of Fukuma-Sakatani [51, 69] it 
is possible to classify the several terms in the entropy current as well as to obtain a first principle 
computation of the total entropy. A related problem is to generalize the action (2.26) to the case where 
the fluid carries either aq = 0, q = loiq = p-brane charge as in the cases studied in [54, 55]. This 
would allow to predict the structure of the piezoeletric moduli measured for charged black branes [6]. 
This problem is related to the entropy current formulation because in both cases it requires obtaining 
from the action a conserved current without additional corrections to the action (2.26) itself. This issue 
is now under investigation and it will be published elsewhere [60]. 

Constraints on the response coefficients: During the analysis of the mode coupling in Sec. 2.6 
we only looked at relations that arise from geometric constraints, such as the Gauss-Codazzi equation 
(2.58). However, there may be stability and thermodynamic constraints imposed by elasticity theory 
or the entropy current analysis in the spirit of [70] that further restrict the set of response coefficients 
obtained in Sec. 2. Furthermore, it would be interesting to understand the physical meaning of each of 
the elastic contributions in (2.40) by obtaining, for example, the corrected speed of propagation of elastic 
waves or the corrections to the elasticity tensor of fluid branes introduced in (2.20). The lack of such 
knowledge is unsatisfactory and deserves further study. 

Elastic corrections to black holes: In Sec. 4 we applied the 2nd order effective action to the case 
of black rings in asymptotically flat space and found that the black hole charges were consistent with 
the free energy interpretation and the Smarr relation. This is compelling evidence that the effective 
description of higher dimensional black holes in the blackfold regime is given by an action of type (2.26) 
with the identification (4.9). It is also important to refer that the predictions made in [5, 6] for the 
corrected horizon angular velocity of (charged) black rings were not accurate enough as they did not take 
into account the contributions to the monopole stress-energy tensor T"^ given in the table of Sec. 2.4. 
We note that we have not obtained a complete prediction of all the corrected black hole thermodynamic 
quantities. This is because we have not developed a first principle computation of the total entropy 
of the system neither accounted for a proper definition of the global temperature in terms of a local 
temperature. Once this is done, the formalism presented here allows for the full predictability of black 
hole charges to order O (e^) in the derivative expansion. We note, however, that due to the lack of 2nd 
order data, at the moment it is only possible to accurately obtain the corrected horizon angular velocity 
for black holes made of bent strings such as the black rings of Sec. 4 and the helical rings and strings 
found in [65]. This is because for embeddings with non- vanishing TZabcd , 2nd order information about the 
response coefficients of the hydrodynamic modes (2.55) is required. In fact, a simple exercise tells us that 
applying the same prescription as for the case of black rings in Sec. 4 to the case of black odd-spheres 
[65] does not lead to consistent results. The study of the bending of the black branes (4.1) to 2nd order 
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would be interesting since it would uncover these extra response coefficients. The bending of strings to 

(e^) is an easier task since the hydrodynamic modes (2.55) will vanish but it would still be a worthy 
endeavour as it would allow to confirm the corrections to the monopole stress-energy tensor found in 
Sec. 2.4 as well as the prediction (4.14). 

We end this work by noting some interesting facts about fluid membranes. When applying the con- 
siderations of Sec. 2 to fluid membranes one does not have to be concerned with backreaction effects 
and hence the effective description of the system is that given by (2.26). If we take the simplest case 
of a membrane embedded in a flat background and only slightly bend it such that TZabcd vanishes, then 
there is a total of 8 response coefficients, 6 of which are related to stationary flows and have not been 
considered previously in the literature. While some of these may be subject to local constraints for real 
cell membranes, as the scalar associated with i?i(k) is, according to the analysis of Sec. 2, this is how 
fluids bend. 
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A Notation, geometry and variations 

In this section we write down in detail the notation used in this work and their relation to the work of 
other authors. We consider surfaces of (p+l)-dimensional worldvolume Wp+i embedded in a background 
space-time endowed with metric 3^2^(3;"). The coordinates x", a = 0, ..., D — \ are space-time coordinates. 
To the worldvolume we assign a set of coordinates a'^, a = 0, ...,p , where p is the total number of spatial 
dimensions of the worldvolume. The worldvolume Wp+i is located at the space-time surface x^ = X^{a"'), 
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where X^(a'^) is a set of mapping functions describing the position of the surface in the ambient space- 
time. Introducing the projector along worldvolume directions Ua'^ = daX^^ we can construct the induced 
metric on the worldvolume as ^ab = g^iuUa^Ub'^, as well as its extrinsic curvature Kab^ = ^ aUb^ symmetric 
in the indices a, h. The covariant derivative along worldvolume directions V^ = u^a^p, compatible with 
both the worldvolume metric jab arid the space-time metric g^^, acts on an arbitrary tensor V'^'^ as 

^aV''' = daV''' + -fab''V'^ + K^uW^ , (A.l) 

where the Christoffel symbols 'yaJ' are computed with respect to the induced metric 7a{, and the Christoffel 
symbols T^^ with respect to the space-time metric gfj_,^. Any space-time vector v'^ can be projected along 
the worldvolume directions using Ua^ such that u" = u"'fj,v^. We further introduce a set of projectors 
n*^, i = 1, ..., D — p — 1 onto the transverse space to the worldvolume defined by g^ui^i^nj'^ = 5ij and 
n^fiUa^ = 0. Any space-time vector can be projected along orthogonal directions such that v^ = n^^v'^. 
The extrinsic curvature by definition is transverse in its third index such that Kab^ = KabUiP. Given 
the normal projectors we can define the extrinsic twist potential via uja^ = —n^^Van^^, which is anti- 
symmetric in the indices i,j. Given the extrinsic twist potential, we can define the outer curvature 
associated with it through the relation [38] 

^ab'^ = VaUJb'^ - VbUjJ^ + UjJ^UJbk^ - UJb'^UJak^ ■ (A.2) 

It is also useful to work with space-time indices and still keeping track of tangential and perpendicular 
components to the worldvolume. For that matter, we introduce the first fundamental tensor ^^^^ such that 
rypy — y^p ^y^iy ^^ab ^ Using 7^^^, any space-time tensor can be projected along the worldvolume while still 
keeping space-time indices. Similarly, we introduce a perpendicular projector such that -L^v= g^u — l^u-, 
which satisfies 7^'^_L^'' = 0. The worldvolume projector Ua^ and the transverse projector n*^ are by 
definition parallel and perpendicular respectively in their space-time index, i.e., Ua^ = Ua'^'j^'^ and n*^ = 
nV-L^/x- Note that both space-time projectors can be defined in terms of the projectors u^a and n'^i such 
that 7^ J, = u^aUv"" and _L^jy = n'^iUi^^. 

Variations d la Capovilla-Guven 

In [38], Capovilla-Guven introduced a covariant derivative V^ that preserves covariance under rotations 
of the normal vectors, such that'^^ 

VaV''=VaV''+u;a'jV''^ , (A.3) 

where Pa is the worldvolume covariant derivative compatible with ^^b- With the definition (A.3), one 
can easily act with it on the transverse indices i but at the expense of introducing torsion. Note that the 
worldvolume covariant derivative Va can be replaced by the covariant derivative Vq introduced in (A.l) 



^"Note that here we are using the opposite conventions compared to the ones used in [38]. That is, the extrinsic curvature 
and extrinsic twist potential defined in [38] should be multiplied by a minus sign to match the conventions used here. 
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if one remembers that Vq does not act on the transverse indices i. In the work presented above we have 
avoided using this terminology by noting that 

VaV' = n\Va {V'nf) = VaV^' - V'^n/Van\ = V^F^^ + Ua\V'''' , (A.4) 

where in the last equality we have used the definition of the extrinsic twist potential. With the definition 
(A. 3) the variation of the extrinsic curvature (2.18) can actually be written as 



S±Kab' = VaVfe$* - R%aj^' - Kac'K%j^^ . (A.5) 

Similarly, for the extrinsic twist potential we have that [38] 

5^u:a'' = -K,b'V'<^'' + KJV'''^'' + R'^ ka^'' . (A.6) 

For variations along the worldvolume directions, the different fields transform with the Lie derivative 
such that covariance under normal rotations is preserved. For example, the tangential variation of the 
twist is 

5iiu;J' = ^''Vbu;J'+u;b''Va^'' , (A.7) 

which when using (A.4) leads to (2.29). An interesting application of this covariant derivative is in the 
rewriting of the spin conservation equation (3.24) which now takes the form 

Vaf'=0 . (A.8) 

We have also considered the variation of the worldvolume Christoffel symbols "/ac as well as of the 
worldvolume Riemann tensor TZabcd- These can be analyzed through the variation of the Christoffel 
symbols, 

Slab"" = ^7"^ i^bhad + Va<57M " V d6jab) ■ (A.9) 

Using this, one has that 

J7^^ed = Ve<57b/ - Vrf^Tfe/ , dnab = Vcd-fab'-Vb6-tac' , (A.IO) 

and finally the variation of the worldvolume Ricci scalar 

5n = Va(l'"'5-fbc^)-V\5-fabn+nab5^''' . (A.ll) 



We have also considered variations of the worldvolume derivative of the extrinsic curvature for co- 
dimension- 1 surfaces. This can be obtained using the chain rule, 



6VaKbc = VadKbc + Sjab''Kdc + djacKbd ■ (A.12) 
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The third fundamental tensor 

In [59] Carter defined the thhd fundamental tensor with space-time indices as 

S.Am' = TA^M^^r'V^/^pa" ' (A-13) 

where the connection V. denotes the covariant derivative along the worldvolunie written with space-time 
indices. We have avoided this notation in the work above by using instead the definition 

V^^^ = 7\VAy^ . (A.14) 

From the third fundamental tensor (A. 13) we can obtain the Codazzi-Mainardi equation simply by taking 
the appropriate anti-symmetrization [59] 

2H[,,]/ = 7^«7^7^^'^ai?p.% . (A.15) 

Contracting this with u^ au'^ b'u^ cn^ u yields Eq. (2.32). 

Variations d la Carter 

For the geodynamic-type models constructed by Carter [36, 37], a difi^erent type of variational principle 
was used than the one used by Capovilla-Guven [38] . It is a Lagrangian variation in which the background 
metric g'^j^ is displaced by an infinitesimal vector such that SlQ^v = 2V(^$,^). We now show that this 
type of variations yields the same results as the ones presented in this work. In order to do so we require 
the Lagrangian variations of the first fundamental tensor 6lJ^'^ = —2ja-^f^'V'^'^'^ and of the second 
fundamental tensor [37] , 



6lK^,P =±^ (V(^V,)Cl>^ - 7^;.7^)^^rp^^ - K\^,)V^^^^ 

Let us consider the simple case of the elastic contribution Ai(k) given in (2.40). The contribution to the 
action is of the form 

nX''] = I V^Xi{k)KfKp , (A. 17) 

where 7 here should be understood as the determinant of 7^jy and k = | — 7^iyk^'k^|-'^'^. We define the 
stress-energy tensor T^'^ and the dipole moment 'D^'^p as before: 

V=7 ^Itiu V=7 oKi_,^P 

where C = ^/—J'Xl(k)KPKp. The Lagrangian variations of (A. 17) are thus of the form 

dL^xn = I V^ (-^T^^'SLif^u + Xi{h)K>'K''6Lg^, + V'^dLKpA . (A.19) 

JWp+i V ^ / 
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The difference here in comparison with (2.30) is the appearance of the middle term above. However, it 
is always cancelled by the last of the terms appearing from the contraction V'^'^ pSiK^^'^ . The end result 
of the variation for all corrections quadratic in the extrinsic curvature is thus^^ 






(A.20) 



The advantage of this approach is that it yields in the end the equations of motion and boundary 
conditions in the form (3.8) and makes it easy to identify the linear momentum V^ . Careful inspection 
of this and comparison with (2.30) leads one to conclude that the variations are the same. For the 
particular case of (A. 17) we have 

Tt"" = -Xi{k)KPKpj>"' + X[{k)ku>'u''KPKp + 4Xi(k)K''^PKp , V^^p = 2Xi(k) j'^^RP , (A.21) 

which when contracted with u"'p,u\, yields the contributions rf and Vf stated in the table given in 
Sec. 2.4. We have mentioned that the middle term in (A. 19) is always cancelled by a contraction involving 
the last term. We briefly show how this is the case. Take the case of the most general action quadratic 
in the extrinsic curvature 

T[X^^]= f ./^\y^-^PK^,"Kxp, , (A.22) 

where the Young modulus yf^^'^P only depends on k, k*" and 7^^". Note that VP^'p = yp'^^'^Kxa^- Now 
introduce a dummie (meaningless in this case) tensor B^''^ defined as 

6^^'^ = ^-^ , (A.23) 

where C = {1/2)^— ^'y'^'^^PRp^'^Kxpa- Then for an action of the type (A.22) this is simply given by 
B^"" = {l/2)y'"^^PK„,,^'Kxp''. Further note that the last term in the contraction T>^"'p6LKf,^P is always 
—VP^PKpy"Va^p- Therefore, BP'^dig^v always cancels this term. 

Linear momentum from the action 

As mentioned in Sec. 3.3 it is possible to obtain the conserved surface currents directly from (A. 17) by 
requiring the action to be invariant under space-time translations along Killing directions. To see this 
note that for an infinitesimal shift along a Killing direction ^P = /3k^ for an arbitrary constant (3. Then 
(A.20) yields 

5lHX^] =13 [ V^ (-V, (V^kx - P^'^^V^kA) + kA^ ^(X^)) . (A.24) 



Note that Carter in [36, 37] has defined the stress-energy tensor T'"' with the opposite sign. 
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Here we have used the definition of Hnear momentum apphed to the case j"'^'" = as given in (3.15). 
Furthermore, we have written the term inside the parenthesis in the second hne of (A. 20) as if^(X'^), 
which when equated to zero yields the equations of motion (3.9)-(3.10) when j"'^'^ = 0. For (A. 24) to be 
invariant under these translations when the equations of motion are satisfied £^{X'^) = we must have 
that 

V, (V^kx - Vf'^V^kx) = . (A.25) 



Indeed, this is recognized as the conservation of the surface current introduced in (3.32) for the case 
jafiu _ Q This is the reason why V^ is called the linear momentum. Moreover, from (A.25), since 
pl'^^] = V.P^t^^l, we have that 






y^p^A _ p^W^A^^^ = ■ (A-26) 



Therefore, along any Killing direction we obtain the equations of motion as written in (3.8). Finally, 
we note that when £'^{X^^) = 0, integrating the first term in (A. 24) into a boundary term yields the 
definition of conserved surface charge (3.36). This type of reasoning had been applied by Guven et al. 
in [61] for extremal branes in flat space. Here we have generalized it for non-extremal branes in curved 
space. 

Conserved charges 

Here we show how to compute the conserved charges from the stress-energy tensor (3.1). Given a spatial 
slice S of the background space-time, the conserved charges are obtained in the usual way 

\Qk\= f d^'-^x^f^^'^k^n, , (A.27) 

here gg stands for the determinant of the spatial part of the background metric g^^. Note that under the 
assumptions of Sec. 3.3 we have that -^—g' = Ro{x°')^/gr and (P^^a^y—^' = Ro{a"')dTdV(^py Further, we 
choose the gauge a^ = t = t, where t is the time coordinate of the background space-time. Hence we 
write (A.27) as 



IQkl = /^ d^'-'^^V^f^^-'j^ = j^ d^'-'xV^T^^'^ V(x") , (A.28) 



where we have defined the tensor field f^uix"') = kf^n^/RQ{x'^). After integrating by parts we obtain the 
final expression 

IQkl = / dV^p^Ro (S^'''V(X^) + B^'^''VpV(X'^)) + boundary term . (A.29) 

The last term is a boundary term which vanishes. Therefore, we are left with the first two terms written 
in a manifestly invariant way under the 'extra symmetry 2' (3.26) and which yield the result (3.38). 
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B Boundary conditions for hydrodynamic modes 

In this appendix we analyze the boundary conditions for the hydrodynamic corrections presented in 
(2.55) and (2.63). We begin with the hydrodynamic scalars given in (2.6). Due to the variation of the 
Christoffel symbols (A. 9) the hydrodynamic scalars contribute with additional terms to the boundary 
equations derived in (2.31). We denote the extra contribution by B^ in terms of 6^ab given in (2.3) and 
S'jab'^ given in (A. 9), which should be added to the last line of Eq. (2.31). 



Scalar 


Ba 


i;i(k)Vi 


rjaviik) (iV^ (kn^n^^Tfee) +7"''J7bcV^k) - f {ku''u'6-f,dV''vi(k) + vi{k)-f^^6-fb,V' 


'k) 


T;2(k)V2 


VaV2ik) ir^d^bc" - ^'"'Slbc'') + Va {r'S7cdV'v2ik) - j''^5^bcV^V2ik)) 


vs{k)Vs 


Va (^V" (i;3(k)k^k^) 5-fbc - V^ (t;3(k)k"k^) 6-fbc + h'"^^bcVb (i;3(k)k'^k^)) - 
rjavsik) {k\^6-fbc'' - k'^k" 6 jbc") 



The complexity of these contributions is useless. For the hydrodynamic corrections (2.6) to fit the pole- 
dipole boundary conditions (3.11), one must require all the above terms to vanish at the brane boundary. 

Boundary contributions due to the quadrupole moment 

We now turn into the case of the boundary contributions of the terms (2.63) due to the presence of the 
quadrupole moment. Besides the boundary terms obtained from (2.30) there is a contribution to the 
variation of the action (2.30) of the form: 



w, 



p+i 



■■\bac 



cba 



•\bac T^d 



6T= / ^/^ValVfdKbc - VbVf^Vc^' - ^'VbVcVT - 2V'""' K"" ^Kbd ^ 



2Vf'' K'^cKbi^r - 2Vf "K^cRbJ^z + V^f'K'^cV^b'^d) + 2Vf''K'^cV^b^a 
2Vf' KW^d^b) + ^'^Vf'^VdKbc + 2V^r^'^VdKdc 



(B.l) 



Note that we used here the index i but one should remember that for co-dimension- 1 surfaces there is only 
one transverse direction. Note also that the first term in (B.l) introduces boundary terms proportional 
to the background Riemann tensor via (2.18). This is a generic effect of quadrupole corrections to the 
equations of motion in the spirit of Sec. 3. Deriving the pole-quadrupole equations of motion is of intrinsic 
interest and it will be published elsewhere. 
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